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ABSTRACT 


The  study  is  dedicated  to  analysis  of  the  information  content  of  outgoing 
microwave  radiation  measurements  and  the  potential  accuracy  of  atmospheric 
temperature  and  humidity  profile  retrieval  for  different  designs  of  experiment. 

It  has  been  shown  that,  in  principle,  outgoing  microwave  radiance  measurements  in 
the  1-200  GHz  spectral  range  contain  a  considerable  quantity  of  independent 
components  of  the  information  (up  to  30)  about  temperature  and  humidity  profiles. 
Extracting  all  these  information  from  outgoing  radiance  spectra  (that  can  be 
achieved  using  optimal  microwave  radiometer)  makes  possible  to  retrieve  the 
temperature  and  humidity  profiles  with  errors  of  0.5-1.0K  and  5-10%  of  the  H20 
mixing  ratio,  respectively. 
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^INTRODUCTION 


At  present,  satellites  are  an  integral  part  of  the  operational  atmospheric 
observing  system.  Remote  sensing  of  meteorological  parameters  is  a  source  of  the 
information  for  analysis  and  numerical  weather  prediction.  In  spite  of  evident 
progress  in  elaboration  of  microwave  atmospheric  sounding  systems  achieved  during 
past  decades,  these  systems  do  not  meet  the  recent  requirements  to  the  accuracy 
and  vertical  resolution  of  temperature  and  water  vapor  profiles,  remotely 
retrieved. 

The  accuracy  and  vertical  resolution  of  remote  sounding  microwave  methods  depend 
on  many  factors: 

1. Measurement  system  characteristics  -  spectral  region  used,  number  and  positions 
of  the  spectral  intervals,  their  spectral  resolution,  instrumental  noise,  absolute 
calibration  quality  of  outgoing  radiance  measurements. 

2. The  accuracy  of  radiative  transfer  microwave  models  used  for  description  of 
radiation  generation  and  transformation  in  the  "atmosphere  -  underlying  surface" 
system  with  due  regard  for  molecular  absorption,  scattering  in  clouds, 
polarization  of  radiation  in  scattering  and  reflecting  processes,  etc. 

3. The  accuracy  of  used  quantitative  characteristics  of  the  interaction  between 
microwave  radiation  and  atmosphere  and  underlying  surface,  namely  -  the  02  and  H20 
absorption  coefficients,  radiative  characteristics  of  clouds  and  surface,  etc. 

4. The  nature  of  the  retrieval  algorithms,  i.e.  the  methods  of  the  inverse  problem 
solution. 

5.  The  quality  and  volume  of  a  priori  information  about  parameters  to  be  retrieved. 

6. The  state  of  atmosphere. 

7. The  methods  of  processing  the  sounding  product  from  satellites. 

All  above  factors  can  influence  on  a  quality  of  retrieved  data  to  the  different 
degree.  But  undoubtedly,  in  our  opinion,  the  measurement  sounding  system 
characteristics  are  the  most '  effective  factor.  In  this  connection,  it  seems  to  be 
very  important  to  analyze  an  information  content  of  outgoing  microwave  radiance 
measurements,  i.e.  to  determine  a  maximal  quantity  of  information  about 
atmospheric  temperature  and  humidity  profiles,  being  in  radiance  measurements  , 
and  to  estimate  the  maximal  potential  accuracy  of  remote  sensing. 

To  solve  these  problems,  there  has  been  developed  a  new  radiative  code  for 
calculations  of  outgoing  brightness  temperature  and  its  variational  derivatives 
with  respect  to  different  atmospheric  and  surface  parameters  in  the  1-200  GHz 
spectral  range  (MIT-SPbSU  program).  The  program  was  successfully  debugged  and  run 
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on  the  PC. 

By  use  of  this  program,  the  dependencies  of  variational  derivatives  (temperature 
and  humidity  weighting  functions)  on  atmospheric  and  surface  states  have  been 
analyzed.  The  calculations  of  the  covariance  matrices  of  outgoing  microwave 
radiance  for  1500  measurement  frequencies  in  the  1-200  GHz  spectral  range, 
simulated  the  "complete  radiance  transfer  experiment",  and  for  SSM/T,  SSM/T-2 
radiometers  of  DMSP  satellite  sounding  system  have  been  carried  out  for  a  case  of 

the  global  statistics  of  the  atmospheric  temperature  and  humidity  profiles.  The 

quantities  of  independent  information  characteristics  contained  in  radiance 

measurements  at  different  measurement  error  levels  for  the  above  designs  of  the 
experiments  have  been  determined  by  the  calculations  of  the  eigenvalues  for 
covariance  matrix  of  outgoing  microwave  radiance. 

The  information  content  of  outgoing  microwave  radiance  measurements  was  estimated 
for  the  complete  radiance  experiment  and  for  SSM/T,  SSM/T-2  radiometers 
(individually  and  in  combination).  The  potential  accuracies  of  atmospheric 
temperature  and  humidity  profile  retrievals  were  calculated  for  different 

measurement  schemes  and  errors,  a  priori  statistics,  vertical  averaging  in  the 
atmosphere. 

The  examples  of  the  numerical  experiments  on  atmospheric  temperature  and  humidity 
profile  retrievals  are  given  for  all  examined  measurement  schemes.  The  above 
studies  have  been  carried  out  for  the  cloudless  atmosphere  and  two  types  of 
underlying  surface  -  sea  and  land. 
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2.RADIATIVE  TRANSFER  MODEL  OF  CLOUDLESS  ATMOSPHERE 
OVER  UNDERLYING  SURFACE  IN  THE  1-200  GHZ  RANGE 

2.1.  Atmosphere  ^ 

These  studies  of  information  content  of  outgoing  microwave  radiative  measurements 
and  potential  accuracy  of  the  temperature  and  humidity  atmospheric  sounding  are 
based  on  the  integral  form  of  the  radiative  transfer  equation.  In  the  case  of  the 
cloudless  atmosphere,  outgoing  radiation  intensity,  a  radiance,  can  be  expressed 
by  the  following  form: 

Po 

J(f)  =  etBfTJP^po)  +  |  BJT(p)]5PffpP))dp  + 

0 

Po 

(l-e^P,(0,p0)-jBfr(p)]-^^dp  (2.1) 

0 

where  ef  is  a  surface  emissivity,  Bf  is  the  Planck  function,  Pf  is  the 
transmittance  function,  Ts  and  T(p)  are  the  surface  temperature  and  the 
atmospheric  temperature  profile,  respectively,  p0  is  a  pressure  near  Earth 
surface. 

For  the  studied  microwave  spectral  range,  using  Rayleigh-Jeans  approach,  the 
equation  (3.1)  can  be  written  in  terms  of  the  brightness  temperature  of  outgoing 
radiance  TB: 

.  Po 

T„(v)  =  efTs-P,(0,p0)  +  jT(p)-6P|(60pP))dp  + 

0 

Po 

(l-e[)-P|(0,p0)-jT(p)6Py)dp  (2.2) 

0 

The  microwave  atmospheric  absorption  spectrum  (and  emission  spectrum)  is 
characterized  by  the  absorption  due  to  molecular  oxygen,  water  vapor,  cloud  and 
precipitation  (Ulaby  et.al.,1982;  Stepanenko  et.al.,1987;  Hollweg  et.al.,  1994). 
Another  minor  gases  in  the  microwave  domain  such  as  ozone,  nitric  acid, 
hydrochloric  acid  and  chemical  radicals  do  not  influence  seriously  on  the 
propagation  of  the  radiation  in  the  case  of  the  nadir  geometry,  but  they  are 
important  for  slant  geometry  (Waters  and  Wofsy,  1978). 
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Molecular  oxygen  strongly  absorbs  due  to  magnetic  dipole  rotational  transitions 
near  60  and  118  GHz.  At  altitudes  below  about  50  km,  the  pressure  broadening  of 
the  lines  dominates  the  Doppler  effect,  which  becomes  important  in  the  middle 
atmosphere.  At  higher  altitudes  (>50  km),  the  Zeeman  splitting  of  the  oxygen 
becomes  important. 

The  water  vapor  lines  in  the  microwave  domain  are  due  to  electric  dipole 
rotational  transitions  extending  into  sub-millimeter  region. 

The  resonant  transition  center  frequencies  for  some  of  the  oxygen  and  water  vapor 
absorption  lines  are  presented  in  table  2.1. 


Table  2.1 

The  resonant  transition  center  frequencies  for  02  and  H20 


(from  Liebe,  1988) 


Line  number 

Water  vapor 

Oxygen 

Absorption  line  center  frequency  (GHz) 

1 

22.235 

tens  of 

2 

183.31 

lines 

3 

323 

at 

4 

325.15 

50-70  GHz 

5 

380.20 

6 

390 

one  line 

7 

436 

at 

8 

438 

118.75 

9 

442 

' 

448 

Non-resonant  (continuum)  interactions  of  the  dry  air  add  a  small  contribution  at 

high  pressure  in  the  lower  atmosphere  due  to  the  Debye  spectrum  of  oxygen  below  10 

GHz  and  pressure-induced  nitrogen  absorption  above  100  GHz  (Mingelgin,1974; 
Stankevich,  1974;  Liebe,  1985,  1989). 

To  accomplish  the  planned  studies  of  the  information  content  and  sounding 
accuracy,  it  is  necessary  to  have  available  of  a  radiative  transfer  model  for  the 
required  spectral  range.  In  the  general  case,  the  radiative  model  is  due  to 
provide  the  calculations  of  the  following  radiation  characteristics  for  arbitrary 

atmospheric  and  surface  models  in  given  spectral  interval  and  with  defined 

spectral  slit  function: 

a) .Outgoing  radiance,  J^f)  and  brightness  temperature  TB(f); 

b) .Variational  (functional)  derivatives  of  outgoing  radiance  (and  brightness 
temperature)  with  respect  to  the  vertical  profiles  and  parameters  of  the 
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atmospheric  and  surface  state  ^ 


x“  A[f’pl 


and 


6TB(f) 


Here  Af  is  considered  spectral  interval,  f  -  its  reference  (usually,  central) 


frequency,  X;  -the  vectors  describing  atmospheric  and  surface  physical  state,  p  - 
pressure  chosen  as  independent  spatial  variable  (as  well  as  in  the  majority  of  the 
studies  on  discussed  problem). 

Thanks  to  kindness  of  Dr.  R.McClatchey  (Phillips  Laboratory,  Atmospheric  Science 
Division),  we  have  available  the  RADTRAN  and  FASCOD3P  codes  permitting 

calculations  of  J^f)  and  TB(f).  In  the  above  radiative  transfer  codes,  the 
blocks  for  variational  derivative  computing  are  absent.  But  keeping  in  mind  the 
opportunities  of  variational  derivative  calculations  by  numerical  methods,  we  have 


tried  to  use  FASCOD3P  and  RADTRAN  codes  for  the  J^f)  an^  TB(f)  calculations. 

In  the  first  stage,  we  would  have  to  make  sure,  that  our  versions  of  these  codes 

give  the  correct  values  of  J^f)  or  TB(f)  and  their  variations  for  different 
atmospheric  models  in  the  given  spectral  range.  As  "a  standard",  we  used,  the  data 
of  International  comparisons  (Chedin  et  al.,  1988).  We  compared  the  results  of 
calculations  using  FASCOD3P,  RADTRAN  codes  and  the  above  data  (Chedin  et  al., 
1988). 

Calculations  of  the  transmittance  functions  and  brightness  temperatures  by 
FASCOD3P  and  RADTRAN  codes  were  carried  out  for  the  following  conditions: 

a) .Spectral  range  from  0.645  to  6.8  cm'1  with  step  of  0.03  cm'1; 

b) .The  vertical  radiation  transfer,  i.e.  calculations  of  the  transmittance  for 
the  whole  of  the  atmosphere  and  the  brightness  temperature  at  the  upper 
atmospheric  boundary; 

c) . Cloudless  atmosphere; 

d) .Surface  temperature  is  equal  to  atmospheric  one  at  the  lowest  level. 

In  the  FASCOD3P  program,  02  and  H20  absorption  and  Voigt  profile  of  spectral 
absorption  line  up  to  64  half-width  have  been  taken  into  account. 

The  comparisons  were  conducted  for  4  atmospheric  models  (1-middle 
latitudes, summer;  2-middle  latitudes,  winter;  3-subarctic, winter;  4-tropical 
atmosphere)  and  2  values  of  the  emissivity  (e=0.5  and  1.0). 

In  the  Table  2.2,  as  an  example,  the  results  of  calculated  transmittance  (P)  and 
brightness  temperature  (TB)  comparisons  for  "middle  latitude,  summer"  model  and 
different  emissivity  are  given. 

In  the  Table  2.2,  the  differences  SP=PF-PR  and  STB=TF-TR  for  a  number  of 
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frequencies  are  given.  The  subscripts  f  and  r  denote  FASCOD3P  and  RADTRAN  codes, 
respectively. 

It  can  be  seen  from  Table  2.2,  for  some  frequencies  the  differences  between  P  and 
Tb  calculations  by  two  codes  may  amount  to  marked  values,  namely,  6P  reaches  to 
0.167,  and  6TB  may  be  larger  than  50K.  For  illustration,  in  Fig.2.1  a)-c)  the 
spectral  dependencies  of  8P  and  8TB  (for  different'  emissivities)  are  given.  The 
ranges  of  maximal  transmittance  and  brightness  temperature  differences  are  located 
at  about  60,  120  GHz  and  at  some  other  spectral  ranges. 

Table  2.2  \ 

Comparison  of  P  and  TB  calculated  by  FASCOD3P  and  RADTRAN 

codes 


Frequency, GHz 

5P 

5TB(e=l) 

8TB(e=0.5) 

19.34 

-0.06395 

-3.12 

-3.45 

25.63 

-0.04479 

-2.47 

-9.47 

31.93 

-0.01712 

-1.35 

-8.84 

36.42 

-0.01791 

-1.45 

-10.28 

45.42 

-0.06707 

-4.79 

-14.62 

47.22 

-0.07530 

-5.88 

-17.76 

51.71 

-0.16176 

-15.94 

-26.44 

53.51 

-0.07556 

-17.48 

-21.42 

55.31 

-0.00008 

1.45 

1.43 

58.01 

0.00000 

11.48 

11.48 

60.71 

0.00000 

12.79 

12.79 

62.51 

0.00000 

-13.61 

-13.61 

64.31 

-0.00008 

0.31 

0.31 

67.90 

-0.13585 

-20.51 

-31.19 

71.50 

-0.02049 

-9.21 

-40.07 

77.80 

0.05063 

-3.21 

-42.04 

84.09 

0.10499 

0.01 

-46.43 

92.19 

0.13508 

1.71 

-50.03 

99.38 

0.15020 

2.49 

-52.57 

105.68 

0.13423 

1.87 

-51.79 

111.97 

0.08774 

-0.30 

-47.93 

117.37 

-0.12776 

-23.37 

-28.52 

119.17 

0.00341 

0.64 

-0.05 

120.07 

-0.10868 

-23.42 

-28.83 

132.66 

0.13579 

1.67 

-52.76 

143.45 

0.13574 

2.02 

-50.85 

155.14 

0.13283 

2.29 

-46.40 

162.34 

0.11527 

2.25 

-39.18 

170.43 

0.06415 

2.02 

-20.45 

173.13 

0.03771 

1.89 

-10.70 

176.73 

0.00589 

1.31 

-0.35 

181.23 

0.00000 

-0.32 

-0.32 

185.72 

0.00000 

-0.04 

-0.04 

198.31 

0.05553 

2.10 

-14.81 

203.71 

0.12962 

3.12 

-27.78 
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difference  tb  difference  Transmittance  difference 


Fig.2.1. Spectral  behavior  of  the  transmittance  (a)  and  brightness  temperature  (b,c) 
differences  at  surface  emissivities  e. 
b-e=1.0;  c-e=0.5. 


The  above  analysis  has  demonstrated,  that  available  FASCOD3P  and  RADTRAN  versions 
give  the  results  differing  very  considerably  for  a  number  cases.  It  testified  to 
the  necessity  of  additional  comparisons  with  independent  data.  In  Table  2.3,  the 
\  comparisons  between  data  calculated  by  FASCOD3P  and  RADTRAN  codes  and  the  results 
of  the  report  about  International  comparisons  of  different  codes  (Chedin  et  al., 
1988)  are  cited. 

Table  2.3 

Comparison  of  monochromatic  brightness  temperatures 
calculated  by  FASCOD3P,  RADTRAN  codes  and  the  results  of 
different  participants  of  the  ITRA 
(Intercomparison  of  Transmittance  and  Radiance  Algorithms) 


Chann. 

AFGL 

LAP 

CNRS 

MIT 

1 

274.80 

274.92 

264.49 

275.03 

2 

274.21 

274.40 

263.48 

274.29 

3 

272.93 

272.90 

263.15 

272.32 

4 

264.99 

264.68 

225.80 

263.21 

5 

250.15 

250.70 

223.67 

260.60 

6 

238.23 

237.80 

222.40 

237.59 

7 

228.47 

228.03 

222.06 

228.58 

8 

222.4 

221.69 

221.77 

221.2 

221.90 

222.16 

9 

218.2 

218.09 

218.31 

217.9 

223.13 

218.14 

10 

219.7 

220.02 

220.41 

219.5 

223.17 

220.43 

11 

223.7 

224.31 

224.45 

223.7 

222.23 

224.77 

12 

230.5 

231.40 

231.70 

230.3 

223.28 

234.50 

13 

241.1 

242.04 

242.23 

240.8 

223.26 

249.95 

14 

253.4 

254.30 

253.58 

253.7 

223.30 

281.33 

15,16 

275.52 

275.90 

255.33 

275.67 

17 

278.14 

277.64 

271.03 

278.30 

18 

245.0 

244.26 

244.34 

244.0 

248.45 

250.39 

19 

258.3 

257.41 

257.58 

257.3 

254.47 

261.14 

20 

271.00 

270.60 

265.64 

271.97 

The  channels  are  numbered  in  accordance  with  AMSU  (Advanced 
Microwave  Sounding  Unit)  channels. 

Data  in  the  table  comes  from  the  following  participants  of  the  Comparisons: 
AFGL  -  Air  Force  Geophysical  Laboratory,  USA; 

LAP  -  Institute  of  Applied  Physics,  University  of  Bern, 

Switzerland; 

CNRS  -  Laboratoire  De  Meteorologie  Dynamique  Du 
Cent.Nat.de  la  Rech.  Scient.,  France; 

MIT  -  Massachusetts  Institute  of  Technology,  USA 


It  is  necessary  to  note  the  evidently  incorrect  results  of  the  brightness 
temperature  (TB)  calculations  for  a  number  of  frequencies  by  the  version  of 
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RADTRAN  code, 

available  to  us: 

For  example: 

Frequency,  GHz 

tb,  k 

56.9678 

428.76 

57.6122 

457.24 

57.6167 

289.96 

It  is  evident,  that  the  first  two  TB  values  make  no  physical  sense. 

The  above  analysis  of  data  calculated  by  available  FASCOD3P  and  RADTRAN  versions 
demonstrates  an  impossibility  of  their  utilization  (or  our  inability  to  use  them) 

for  our  studies. 

Therefore  we  developed  the  new  own  code  MIT-SPbSU  for  J^f)  and  TB(f) 

calculations  on  the  basis  of  detailed  radiative  transfer  model  MIT  description 

(Liebe  et  al.,  1993).  Simultaneously  we  have  not  only  repeated  the  Liebe 

algorithm,  but  also  have  added  a  block  for  variational  derivative  calculations. 
The  latter  results  in  the  suitability  of  our  program  for  inverse  problem  modeling. 
Radiative  transfer  model  MIT  is  one  of  the  latest  version  of  MPM  code  (Millimeter 
wave  Propagation  Model),  most  widely  used  in  practice  and  originally  described  by 
Liebe  (1981).  It  has  been  updated  several  times  (Liebe,  1985,  1988,  1989;  Liebe 
and  Hufford,  1989a;  Liebe  et  al.,  1989b;  Liebe  et  al.,  1991;  Liebe  et  al., 

1993)  to  match  spectroscopic  measurements  of  absorption,  propagation  delay  and 
emission  of  molecular  oxygen,  water  vapor  and  atmosphere.  Liebe  (MPM)  integrates 
oxygen  and  water-vapor  lines  up  to  1  GHz  assuming  Van  Vleck-Weisskopf  line  shapes 
with  line  overlap  corrections  according  to  Rosenkranz  (1975).  He  adds  dry  air  and 
water-vapor  continue  which  are  mainly  based  on  own  laboratory  measurements.  A 
strong  negative  temperature  dependence  of  the  continuum  absorption  is  included. 
Below  300  GHz  a  modified  continuum  reduces  the  line  parameter  base  to  a  smaller 
amount  of  spectroscopic  parameters  to  be  considered  thus  reducing  the 

computational  effort  when  calculating  transmittances  and  radiances. 

MIT  code  (Liebe  et  al.,  1993)  models  of  the  propagation  characteristics  of  the 
atmosphere  for  the  frequency  range  from  1  to  1000  GHz  (ITHz).  Refractivity  spectra 
of  the  main  natural  absorbers  (i.e.  oxygen,  water-vapor,  suspended  droplets  and 
ice  particles)  are  computed  from  known  meteorological  variables.  The  primary 
contributions  of  dry  air  come  from  44  02  lines.  Results  from  extensive  60  GHz 
laboratory  measurements  of  the  pressure-broadened  02  spectrum  were  applied  to 
update  the  line  data  base.  The  water-vapor  module  considers  34  local  H20  lines 
plus  continuum  contributions  from  the  H20  spectrum  above  1  THz,  which  are 
formulated  as  wing  response  of  a  pseudoline  centered  at  1.8  THz. 
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The  authors  have  maintained  that  the  new  code  MPM-93  reproduces  the  spectral 
characteristics  of  the  clear  atmosphere  (02,  H20)  between  18  and  930  GHz  within 
the  in  certainty  limits  of  five  reported  controlled  experiments  (Becker  and 
Autler,1946;  Liebe,  1987;  Furashov  et  al.,  1989;  Bauer  and  Godon,  1991;  Godon 
et  al.,  1992). 

The  comparisons  of  brightness  temperatures  calculated  by  MIT-SPbSU  program  with 
the  results  represented  in  the  report  (Chedin  et  al.,  1988)  for  a  number 
frequencies  are  given  in  Table  2.4. 

Table  2.4 

Comparisons  of  data  calculated  by  MIT-SPbSU  program  with 
the  results  of  the  report  (Chedin  et  al.,  1988) 


Chann. 

AFGL 

LAP 

CNRS 

MIT 

MIT-SPbSU(l) 

MIT-SPbSU(2) 

1 

274.80 

274.92 

274.69 

274.69 

2 

274.21 

274.40 

274.15 

274.15 

3 

272.93 

272.90 

272.41 

272.41 

4 

264.99 

264.68 

264.31 

264.31 

5 

250.15 

250.70 

249.62 

249.62 

6 

238.23 

237.80 

237.82 

237.82 

7 

228.47 

228.03 

228.36 

228.36 

8 

222.4 

221.69 

221.77 

221.2 

221.74 

221.75 

9 

218.2 

218.09 

218.31 

217.9 

218.30 

218.30 

10 

219.7 

220.02 

220.41 

219.5 

220.33 

220.33 

11 

223.7 

224.31 

224.45 

223.7 

224.77 

224.77 

12 

230.5 

231.40 

231.70 

230.3 

232.14 

232.15 

13 

241.1 

242.04 

242.23 

240.8 

242.83 

242.83 

14 

253.4 

254.30 

253.58 

253.7 

253.94 

253.94 

15,16 

275.52 

275.90 

275.45 

275.43 

17 

278.14 

277.64 

277.18 

277.16 

18 

245.0 

244.26 

244.34 

244.0 

243.89 

243.89 

19 

258.3 

257.41 

257.58 

257.3 

256.51 

256.51 

20 

271.00 

270.60 

269.53 

269.53 

In  the  tab.2.4,  the  monochromatic  value  for  nonmonochromatic  channels  is  the  mean 
of  all  brightness  temperatures  calculated  for  the  sideband  center  frequencies.  MIT 
SPbSU(l)  and  MIT  SPbSU(2)  are  the  values  of  radiance  and  brightness  temperature, 
respectively,  calculated  by  (1.1)  and  (2.2).  Other  designations  are  in  keeping 
with  those  in  tab.2.3. 

It  can  be  seen  that  the  accordance  of  our  results  with  the  independent  data  is 
very  satisfactory.  As  a  rule,  the  differences  do  not  exceed  IK,  and  those  are  the 
tenths  of  Kelvin  in  the  most  cases. 

Developed  MIT-SPbSU  code  was  further  used  for  studying  the  information  content  of 
outgoing  microwave  radiation  measurements  and  to  analyze  a  remote  sensing  accuracy 
under  different  designs  of  satellite  experiments. 
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2.2.Surface  emissivity 

Most  emissivity  models  used  in  practice  are  empirical  or  semiempirical  ones.  The 
review  of  these  model  has  been  given  in  our  "Feasibility  Study"  (Timofeyev  et  al., 
1994). 

Though  the  use  of  theoretical  models  is  more  desirable,  in  the  most  cases  their 

applying  is  impossible  as  their  mathematical  formulations  are  very  complex  and 
actual  values  of  employing  numerous  parameters  are  unknown. 

Pure  empirical  models  are  developed  by  fitting  different  functions  of  measurable 

ground  truth  parameters  into  the  observed  emissivity  data.  Unfortunately,  as  we 
have  mentioned  in  the  "Feasibility  Study"  (Timofeyev  et  al.,  1994),  experimental 
emissivity  data  in  the  examined  frequency  range  are  very  limited.  Semiempirical 
models  use  also  some  theoretical  dependencies  in  addition  to  emissivity 

observations. 

Since  empirical  models  are  received  from  measurements  at  certain  conditions 
(frequencies,  polarization,  etc.),  their  usability  is  restricted.  By  applying 

semiempirical  models,  the  range  of  their  practical  implementation  may  be 
enlargened  to  necessary  measurement  conditions. 

To  estimate  a  range  of  emissivity  variations  for  different  natural  surfaces,  we 

have  used  the  semiempirical  models  of  land  and  sea  in  available  microwave 

frequency  domain  and  have  developed  the  computer  programs  for  calculations.  Only 
summer  conditions  have  been  examined  . 

2.2.1.Sea  emissivity 

For  calm  sea  surface,  a  penetration  depth  is  less  than  a  half  centimeter  at  the 

lowest  examined  frequency  (Ulaby  et  al.,1986),  and  much  less  in  the  higher  ones. 
So  emission  behavior  of  open  sea  water  in  the  absence  of  foam  is  exclusively 
determined  by  the  surface  layer  and  can  be  calculated  from  the  power  reflection 
coefficients  (Fresnel  ones)  of  the  surface.  Their  dependence  on  salinity,  surface 
temperature  and  frequency  have  been  well  defined  (e.g.  Stogryn,1971;  Klein  and 
Swift, 1977). 

The  main  conclusion  from  the  studies  of  calm  sea  surface  emissivity  is  the 

following.  Above  6  GHz,  the  only  effecting  parameter  (in  addition  to  angle  of 
incidence,  polarization  and  operation  frequency)  is  the  surface  temperature. 

The  emissivity  of  wind  roughened  sea  surface  increases,  that  is  caused  by  surface 
roughness  and  foam  cover  appearing  on  surface  at  wind  speed  of  about  7m/c.  The 
foam  cover  influence  on  emissivity  prevails  and  depends  on  wind  speed.  For  foam 


11 


cover  influence  accounting,  we  used  a  parameterization  proposed  by  Wilheit  (1979) 
and  developed  by  Wilheit  and  Chang(1980). 

Wilheit  model  can  be  applied  up  to  37  GHz.  But  since  the  model  is  based  on  the  use 
of  a  real  surface  geometry,  it  probably  can  be  applied  within  the  whole  of 
examined  frequency  range.  Besides,  a  comparison  of  different  ocean  surface 
emissivity  models  with  measurements  given  in  the  papers  of  Sasaki  et  al.  (1987; 
1988)  confirmed  that  Wilheit  model  seems  to  be  more  accurate  at  the  higher 
frequencies  (Pulliainen  et  al.,  1990). 

In  Fig.2.2  the  data  of  sea  surface  emissivity  calculations  are  given  for  calm  sea 
and  sea  driven  by  wind  with  speeds  of  20  and  40  m/c.  It  can  be  seen,  that  the 
variations  of  sea  surface  emissivity  are  fairly  large  and  range  from  0.07  at  3  GHz 
to  0.2  at  43  GHz.  It  seems  possible  to  parameterize  a  spectral  dependence  of  sea 
emissivity  at  high  frequencies  by  linear  function  for  inversion  purposes.  This 
conclusion  coincides  with  the  result  of  formerly  study  (Ferraro  et  al.,  1986). 

2.2.2.Land  emissivity 

To  estimate  the  land  emissivity,  we  have  used  the  models  of  surface  differing  with 
soib  texture  (sand  or  clay),  categories  of  roughness  (smooth;  low,  medium  or  high 
roughness)  and  types  of  vegetation  (short  grass,  agricultures,  forest). 

The  emissivity  of  bare  smooth  soil  is  determined  by  dielectric  properties  of  the 
target.  Spectral  behavior  of  soil  dielectric  constants  depends,  in  the  main,  on 
soil  moisture  and,  to  small  degree,  on  texture. 

For  our  calculations  we  have  used  the  laboratory  data  on  soil  dielectric  constant 
given  by  Leschansky  et  al.  (1971). 

Frequency  dependencies  of  bare  smooth  sand  and  clay  surface  emissivity  at 
different  moistures  are  given  in  Fig.2.3.  As  seen,  the  emissivity  variations, 
depending  on  volumetric  soil  moisture  increase  from  5  to  20%,  range  from  0.05  at 
low  frequencies  to  0.12  at  high  ones. 

Surface  roughness  results  in  an  increase  in  emissivity  over  that  of  an  equivalent 
smooth  surface.  One  of  the  first  practical  quantitative  semiempirical  approaches 
to  dealing  with  surface  roughness  was  presented  Choudhury  et  al.  (1979). 
Emissivity  of  roughness  surface  is  calculated  from  smooth  surface  emissivity  using 
an  empirically  determined  roughness  parameter,  that  is  proportional  to  the  RMS 
height  variations  of  the  surface.  The  variants  of  this  approach  were  proposed  by 
Wang  et  al.  (1980)  and  Konyate  et  al.  (1981).  More  sophisticated  models  have  been 
also  studied  but  are  lacking  in  verification  or  practicality. 

Developing  our  computer  program,  we  have  preferred  the  model  of  Wang  et  al.  (1980) 
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emissivity 


Fig. 2. 2. Spectral  dependence  of  sea  surface 

emissivity  at  different  wind  speeds  v. 
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Fig. 2. 3. Spectral  dependence  of  bare  smooth  land 
emissivity  at  different  volumetric  soil 
moisture  w. 
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because  of  its  better  correspondence  to  experimental  data  (Promes  et  al.,  1988; 
Tosellini,  1989). 

Influence  of  surface  roughness  on  emissivity  of  dry  and  wet  sand  is  given  in 
Fig.2.4.  For  dry  fields,  the  effect  of  surface  roughness  on  emissivity  is  small, 
but  increases  with  soil  moisture  growth.  Emissivity  variations  are  about  0.04  and 
reach  to  0.1  for  dry  and  wet  soils,  respectively.  At  high  frequencies  (above  10 
GHz),  the  spectral  behavior  of  emissivity  can  be  simulated  as  linear  dependence. 

A  vegetation  canopy  over  the  soil  increases  the  emissivity  very  strongly.  A 
theoretical  model  of  this  process  that  treats  the  problem  by  using  the 

electromagnetic  relationships  of  a  two  layer  incoherent  nonscattering  medium  was 
developed  by  Basharinov  and  Shutko  (1975).  In  this  model  the  vegetation  is  treated 
simply  an  absorbing  layer.  The  nature  of  the  attenuation  of  this  layer  was  modeled 
by  Kirdiashev  et  al.  (1979)  as  optical  depth  and  vegetation  wet  biomass.  The  plant 
shape  parameter  is  empirical.  Naturally,  there  are  some  other  studies  of 

vegetation  canopy  emissivity  (for  example,  Jackson  and  Schmugge,  1991;  Jackson  and 
O’Neill,  1990),  but  only  the  study  of  Kirdiashev  give  all  the  data  needed  to 
develop  the  software. 

Having  in  mind  that  the  natural  surfaces  are  usually  characterized  by  both 
vegetation  canopy  and  roughness,  we  have  developed  also  a  joint  model. 

In  Fig.2.5,  as  an  example,  we  give  the  calculated  spectral  behaviors  of  surface 
emissivity  for  dry  and  wet  roughness  land  covered  with  vegetation  canopy  (short 
grass).  For  comparison,  the  spectral  dependencies  of  bare  dry  and  wet  soil 
emissivity  are  given  in  Fig.2.5  too.  It  can  be  seen  from  Fig.2.5,  that,  as  it  was 
in  the  previous  cases,  emissivity  can  be  presented  by  linear  spectral  dependence 
at  frequencies  exceeding  10  GHz. 

Summing  the  above  estimates  of  emissivity  variability  range,  it  can  be  noted  the 
following: 

1. Emissivity  models  of  different  surfaces  have  been  developed,  tested  and  work 
well  only  up  to  about  40GHz.  At  higher  frequencies,  the  influence  of  scattering 
processes  on  emissivity  of  both  sea  and  roughness  and  vegetated  land  is  very 
expected.  Thus,  further  theoretical  and  experimental  studies  are  needed  to  enhance 
the  models  at  higher  frequencies. 

2. The  case  of  calm  sea  surface  is  the  most  clear.  The  frequency  (f)  dependence  of 
calm  sea  emissivity,  e,  can  be  given  by  linear  expression  (Chedin  et  al.,  1988): 

e  =  0.362  +  0.00272  f[GHz]  (2.3) 

3. Emissivity  variations  caused  by  wind  influence  can  reach  up  to  0.2  at  f=40  GHz 
and  decrease  with  frequency  growth. 
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emissivity 


Fig. 2. 4. Effect  of  surface  roughness  on  emissivity 
of  dry  (w=5% )  and  wet  (w=20& )  sand. 


16 


emissivity 


Fig. 2. 5. Effect  of  vegetation  canopy  (short  grass)  and 

medium  roughness  on  emissivity  of  dry  (w=5%) 
and  wet  (w=20%)  sand. 
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4. Land  eraissivity  ranges  from  0.7  to  0.95  depending  on  soil  moisture,  roughness 
and  vegetation  cover. 

5. At  frequencies  exceeding  10  GHz,  spectral  behavior  of  different  underlying 
surfaces  can  be  simulated  by  linear  expressions  with  different  slopes. 

It  is  naturally,  that  for  atmospheric  inversion  purposes,  we  ought  not  model  these 
surface  variations  in  frill  measure.  In  our  inversion  study  as  a  first  step,  we 
simulate  a  spectral  sea  emissivity  dependence  by  expression  (2.3)  to  have  a 
possibility  of  comparing  our  retrieval  study  with  previous  ones.  Land  emissivity 
is  supposed  to  be  constant  and  equal  to  0.95.  The  choice  of  this  large  value  of 
the  emissivity  is  conditioned  by  desire  of  well  distinguishing  of  sea  and  land. 

23.Examples  of  variational  derivative  calculations  by  MIT-SPbSU  code 


To  demonstrate  the  potential  opportunities  of  developed  MIT-SPbSU  code  and  to 
analyze  the  connections  between  outgoing  brightness  temperature  TB(f)  and 
atmospheric  temperature  T(z)  or  humidity  q(z)  profiles,  we  give  the  examples  of 
calculations  of  TB(f)  variational  derivatives  with  respect  to  T(z)  and  q(z). 

Variational  derivatives  characterize  (in  linear  approach)  the  sensitivity  of 
outgoing  radiance  (the  brightness  temperature)  to  atmospheric  parameter 
(temperature,  humidity)  variations  at  different  atmospheric  levels.  It  follows 
from  Taylor-series  expansion  of  STB(f): 


Po 


8TB(f)  =  j 


6T„(f) 

5T® 


Po 


T(p) 


r  6TB(f) 

■6T®de  +  J  Sqfj 


0  0 
which  can  be  rewrite  in  the  operator  form  as: 

6TB(f)  =  At6T  +  AqSq 


q(p) 


■6q(p)dp 


(2.3) 


(2.4) 


It  is  necessary  to  note,  that  variational  derivatives  differ  from  so-called 
"weighting  functions"  used  in  some  works  concerning  microwave  spectral  range  (e.g. 
(Isaacs  et  al.,  1987).  The  correct  expressions  for  these  variational  derivatives 
has  been  given  in  paper  of  Timofeyev  (1974).  Similar  "weighting  functions"  (those 
would  be  more  properly  termed  as  the  variational  ■'derivatives)  were  applied  by  Kuo 
et  al.(1994). 

Because  of  numerous  wrong  representations  especially  of  humidity  weighting 
functions,  we  adduce  here  the  expression  for  variational  derivative  with  respect 
to  atmospheric  humidity  q(z)  (Timofeyev,  1974). 

In  linear  approach,  the  variations  of  TB  (6TB),  conditioned  by  the  deviations  of 
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q(z),  6q,  from  some  average  q,  may  be  expressed  as: 
Po 

6TJf)  I 

•6q(p)dp 


=  J  SqfpJ 


where 


6TE(f) 

«q(p) 


0 


q(p) 


q(p) 

=  Aq  is  humidity  variational  derivative. 


(2.5) 


By  use  of  (2.3),  one  can  receive  for  A 


\  IAq> 


i=l 


A„(P)  =  -  e(T,  -  T(Pa)} 


P(Po)k(p) 

g 


Po 


A,2(P)  =  -  P(P')dp' 

P 

P2(Po) 


Aq3(p)  =  2(l-e)T(p0)k(p) 


g 


AqjCp)  =  (l-e)k(p)^J^M  P*(p)dp 

0 

A,5(P)  =  (l-e)k(p)?|^J  ^  P*(p)dp', 


where 


P(p)  =  exp 


P  (p)  =  exp 


r 

|  J  k(p)q(p)dp 


if 


k(p)q(p)dp 


(2.6) 


k(p)  -  the  absorption  coefficient,  e  -  the  emissivity,  Ts  -  the  surface 
temperature,  g  -  free  fall  acceleration,  q(p)  -  specific  atmospheric  humidity. 

In  Fig.2.6,  there  are  shown  the  monochromatic  TB(f)  variational  derivatives  with 
respect  to  T(z)  for  7  frequencies  corresponding  to  SSM-T  channels.  Fig.2.6 
demonstrates  the  vertical  opportunities  of  temperature  sounding  using  outgoing 
radiative  measurements  in  the  02  band  and  fairly  uniform  covering  of  all  height 
ranges  by  variational  derivatives  for  different  frequencies. 
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dTB/dT(p)  [K/(K*mb)] 


Fig. 2. 6. Variational  derivatives  of  brightness  temperature 
with  respect  to  atmospheric  temperature  profile. 
Atmospheric  model  -  middle  latitudes,  summer; 
surface  emissivity  is  0.95  (land); 

Ts=Tioi3=294.2K. 
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It  is  worth  noting  a  growth  of  maximal  derivative  values  with  pressure  decrease. 
It  is  caused  by  choice  of  pressure  as  a  parameter.  To  remove  an  illusory 
unequivalence  of  the  variational  derivatives,  those  are  to  be  multiplied  by  the 
pressure.  Then  the  variational  derivatives  will  be  illustrate  the  influence  of 
parameter  variations  in  the  layers  of  equal  thickness  (regardless  of  pressure)  on 

brightness  temperature. 

To  illustrate  the  above  arguments  let  us  compare  Fig.2.6  with  Fig.2.7  representing 
the  same  derivatives  but  multiplied  by  corresponding  pressures.  It  is  evident, 
that  the  latter  figure  is  more  informative  and  all  following  data  concerning  both 
temperature  and  humidity  profiles  will  be  expressed  in  terms  of  this  product, 

named  by  proper  "weighting  functions". 

To  illustrate  surface  influence  on  TB(f)  derivatives,  Fig.2.8  gives  a  comparison 

the  temperature  weighting  functions  calculated  with  the  values  of  surface 
emissivities  typical  for  land  and  sea.  The  differences  have  been  observed  only  for 
three  examined  frequencies,  therefore  in  Fig.2.8  the  rest  curves  are  absent.  As  it 

was  expected,  the  surface  influence  on  the  temperature  weighting  functions  is  the 

larger  than  the  derivative  maximum  is  the  nearer  to  the  Earth  surface.  The  curve 

differences  are  observed  only  for  the  pressures  which  are  larger  than  100  mb,  that 
let  us  conclude  that  the  uncertainties  in  surface  emissivity  will  not  influence  an 

accuracy  of  temperature  retrieval  in  higher  atmospheric  layers. 

It  is  paid  attention  an  extraordinary,  at  first  sight,  behavior  of  the  temperature 

weighting  functions  with  negative  values  near  surface  for  the  frequency  of  50.5 
GHz.  This  appearance  is  caused  by  a  strong  temperature  dependence  of  water  vapor 
absorption  near  Earth  surface.  In  this  spectral  range,  the  absorption,  and 

atmospheric  radiance  respectively,  decreases  abruptly  with  temperature  increase. 

When  surface  emissivity  is  equal  to  0.95,  this  effect  is  masked  by  dominant 

contribution  of  surface  emission  into  outgoing  radiance  and,  as  a  consequence, 

into  temperature  weighting  function.  It  is  worth  to  note,  that  such  strong 

dependencies  of  variational  derivatives  on  the  atmospheric  state  testify  to  poor 

applicability  of  purely  linear  approach  for  the  solution  of  the  inverse  problem 
and  to  the  necessity  of  using  the  iterations. 

Fig.2.9  describes  an  influence  of  the  surface  temperature  or  the  difference 

between  surface  and  atmospheric  near  surface  layer  temperatures  (AT)  on 

temperature  weighting  functions.  The  results  are  given  for  two  frequencies  only, 

as  the  influence  is  not  observed  for  other  examined  ones.  It  is  seen  the 

negligible  variation  of  derivative  values  even  for  temperature  differences  of  10K. 

To  analyze  the  influence  of  different  atmospheric  parameters  on  the  temperature 
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Pressure  [mb] 


Fig. 2. 7. Temperature  weighting  functions  for 
different  frequencies. 

Atmospheric  model  -  middle  altitudes,  summer; 
surface  emissivity  is  0.95  (land); 

T8=T1013=294K. 
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Fig. 2. 8. Influence  of  surface  emissivity  on  temperature 
weighting  functions  for  different  frequencies. 
Atmospheric  model  -  middle  latitudes,  summer 
surface  emissivity:  0.95  -  land,  0.5  -  sea; 
T„=T,m,=294.2K. 
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Fig. 2. 9. Influence  of  surface  temperature  on  temperature 
weighting  functions  for  different  frequencies. 
Atmospheric  model  -  middle  latitudes,  summer; 
surface  emissivity  -  0.95  (  land); 

Tioi3=294.2K. 
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weighting  functions,  we  carried  out  the  calculations  for  a  number  of  simulated 
T(z)  and  q(z)  sets.  Fig.2.10  illustrates  the  influence  on  temperature  weighting 
functions  caused  by  changing  the  summer  T(z)  for  middle  latitudes  by  winter  one, 
Fig.2.11  and  2.12  show  the  effect  of  the  same  change  of  q(z),  and  both  T(z)  and 
q(z)  profiles,  respectively. 

Joint  analysis  of  these  three  figures  makes  possible  to  conclude  the  following: 

1. Temperature  profile  variations  slightly  influence  the  temperature  weighting 

functions  because  of  weak  atmospheric  transmittance  dependence  on  the  temperature. 

2. Humidity  profile  variations  effect  the  temperature  weighting  functions  in  the 
layers  below  pressure  level  of  100  mb,  and  relative  variation  of  derivative  values 
may  achieve  to  about  20%. 

3. Under  joint  changes  of  temperature  and  humidity  profiles,  these  effects  are 
added  together. 

As  a  whole,  the  fulfilled  analysis  let  us  conclude  that  temperature  weighting 
functions  slightly  depend  on  different  parameters  of  the  atmospheric  layers  with 
the  exception  of  the  near  surface  atmospheric  layers  characterizing  by  the  strong 
influence  of  surface  emissivity  and  atmospheric  humidity  on  microwave  radiance 
transfer.  It  means,  that  the  methods  based  on  linearization  may  be  used  (taking 
into  account  the  above  remark  about  near  surface  layers)  for  atmospheric 
temperature  profile  retrieval. 

Analyzing  the  humidity  variational  derivatives  (variational  derivatives  with 
respect  to  humidity  profile)  we  dealt  with  those  multiplied  by  the  pressure  values 
as  it  was  justified  previously.  In  addition,  it  is  more  visual  to  deal  with 
relative  water  vapor  variations  at  each  level  instead  of  with  absolute  atmospheric 
water  vapor  content.  This  approach  is  justified  by  the  proportionality  of  water 
vapor  variability  to  average  value  of  mixing  ratio  (which  is  strong  varies  with 
altitude)  at  each  level.  In  Fig.2.13  and  2.14  to  illustrate  these  aspects,  we  give 
the  variational  derivatives  of  outgoing  brightness  temperature  with  respect  to 
water  vapor  mixing  ratio  profile  and  those  multiplied  by  the  pressure  and  H20 
mixing  ratio  and  named  as  "relative  humidity  weighting  functions".  The  weighting 
functions,  represented  here  and  later,  are  calculated  for  frequencies 
corresponding  to  three  SSM/T-2  channels  located  near  the  H20  absorption  line  of 
183.310091  GHz,  to  two  its  channels  near  90  and  150  GHz  and  to  the  center  of  the 
above  line.  Let  us  analyze  an  influence  of  different  parameters  on  variational 
derivatives  of  outgoing  brightness  temperature  with  respect  to  H20  content 
profile. 

The  influence  of  surface  emissivity  appears  only  for  two  examined  frequencies, 
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Pressure  [mb] 


Fig. 2. 10. Influence  of  atmospheric  temperature  profile 

on  temperature  weighting  functions  for  different 
frequencies. 

Atmospheric  model  -  middle  latitudes,  summer; 
surface  emissivity  -  0.95  (land); 

Ts  =  T1013. 

Atmospheric  temperature  profile  corresponds  to 
summer  or  winter  model  for  middle  latitudes. 


26 


Pressure  [mb] 


dTB/dT(p)*p  [K/K] 


Fig. 2. 1 1. Influence  of  atmospheric  humidity  profile  on  temperati 
weighting  functions  for  different  frequencies. 
Atmospheric  model  -  middle  latitudes,  summer; 
surface  emissivity  -  0.95  (land); 

Ts^Tjoia. 

Atmospheric  humidity  profile  corresponds  to 
summer  or  winter  model  for  middle  latitudes. 
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Fig.2.13.Variational  derivatives  of  brightness  temperature 

with  respect  to  water  vapor  mixing  ratio  variations 
for  different  frequencies. 

Atmospheric  model  -  middle  latitudes,  summer; 
Ts=T10i3=294K. 
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Fig. 2. 14. Relative  humidity  weighting  functions  for 
different  frequencies. 

Atmospheric  model  -  middle  altitudes,  summer; 
Ts=T10i3=294K. 


corresponding  the  relative  humidity  weighting  functions  are  given  in  Fig.2.15.  It 
may  be  seen  the  strong  dependence  of  the  relative  humidity  weighting  functions  on 
the  surface  emissivity.  Their  values  increase  in  three,  when  surface  emissivity 
varies  from  0.95  to  0.75.  This  effect  is  caused,  evidently,  by  increasing  the 
atmospheric  contribution  into  outgoing  radiance  with  decreasing  the  surface 
emissivity  and  may  be  analyzed  quantitatively  on  the  basis  of  Eqs.  (2.5)-(2.6). 

Fig.2.16  demonstrates  also  a  marked  dependence  of  the  relative  humidity  weighting 
functions  on  the  difference  between  near  surface  atmospheric  temperature  and 
surface  one  (AT).  The  increase  of  this  difference  from  5  to  10K  results  in 
relative  growth  of  derivative  values  by  tens  of  percentages. 

To  analyze  an  influence  of  different  atmospheric  parameters  on  the  relative 
humidity  weighting  functions,  we  have  fulfilled  the  calculations  for  various  sets 
of  simulated  T(z)  and  q(z)  profiles. 

Fig.2.17  illustrates  the  influence  on  the  relative  humidity  weighting  functions 
due  to  a  change  of  summer  T(z)  for  middle  latitudes  by  winter  one,  Fig.2.18  and 
2.19  show  the  effect  of  the  same  change  of  q(z),  and  both  T(z)  and  q(z)  profiles, 
respectively. 

Joint  analysis  of  these  three  figures  makes  possible  to  conclude  the  following: 

1. Temperature  profile  effects  the  relative  humidity  weighting  functions  noticeably 
for  all  three  examined  channels.  At  some  levels,  the  derivative  values  may  vary  by 
the  factor  of  two. 

2. Humidity  profile  influences  the  relative  humidity  weighting  functions  very 
strongly.  It  is  observed  the  variations  of  derivative  values  by  several  orders  of 
magnitude,  the  shifts  of  derivative  maxima  by  values  about  100  mb  and  the  change 
of  derivative  signs  due  to  humidity  profile  variations. 

3. The  joint  variations  of  T(z)  and  g(z)  amplify  mutually  these  effects. 

These  results  demonstrate  that  inverse  problem  with  respect  to  humidity  profile  is 
strongly  nonlinear. 
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Fig. 2. 16. Influence  of  surface  temperature  variations  (  T) 
on  relative  humidity  weighting  functions  for 
different  frequencies. 

Atmospheric  model  -  middle  latitudes,  summer; 
Tioi3=:294K. 
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Fig. 2. 17. Influence  of  atmospheric  temperature  profile 
on  relative  humidity  weighting  functions  for 
different  frequencies. 

Atmospheric  model  -  middle  latitudes,  summer; 

ts=t1013. 

Temperature  profile  corresponds  to  winter 
or  summer  model  for  middle  latitudes. 
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Fig. 2. 18. Influence  of  atmospheric  humidity  profile  on  relative 
humidity  weighting  functions  for  different  frequencies. 
Atmospheric  model  -  middle  latitudes,  summer; 

Ts=T10l3. 

Humidity  profile  corresponds  to  summer  or  winter 
models  for  middle  latitudes. 
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Fig. 2. 19. Influence  of  atmospheric  model  on  relative  humidity 
weighting  functions  for  different  frequencies. 

Atmospheric  models-middle  latitudes,  summer  or  winter; 

ts=t1013. 
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3.CALCULATIONS  OF  INFORMATION  CONTENT  OF  OUTGOING 
MICROWAVE  RADIANCE  MEASUREMENTS 

3.1.Information  content  of  remote  sensing  measurements 

To  compare  the  different  spectral  designs  of  the  measurements,  it  is  necessary  to 

know  the  basics  of  the  quantitative  methods  for  estimation  and  comparison  of  the 
measurement  efficiencies. 

Let  us  consider  the  most-used  lines  of  relevant  approaches. 

The  inverse  problem  of  remote  atmospheric  sensing  is  a  solution  of  the  operator 

equation: 

y  =  A-x  +  £  (3.1) 

where  x  is  the  vector  of  sought  parameters  (atmospheric  temperature  and  humidity 
profiles  and  surface  temperature  in  the  case  discussed);  y  is  the  vector  of 

quantities  to  be  measured  by  spectral  device;  e  is  the  vector  of  measurement 
noise.  The  operator  A  describes  the  radiative  transfer  process  and  links  the 

atmospheric  state  parameters  Xj  with  radiative  characteristics  yj  measured  by  a 
device.  (i=l,2,...m;  j=1.2...n;  where  m  is  the  number  of  parameters,  describing 
the  atmospheric  stale;  n  is  the  number  of  the  measurements.) 

In  this  treatment  of  the  linear  statistical  method,  A  is  assumed  to  be  linear. 

The  every  sets  of  the  sought  atmospheric  parameters  will  be  considered  to  belong 
to  the  vector  space  X.  And  the  vector  space  Y  will  be  considered  to  contain  the 
sets  of  the  measured  values  and  measurement  noise. 

The  noise  e  is  usually  assumed  to  be  distributed  normally  with  the  ensemble 
average  i  and  covariance  matrix  2: 

e  =  Ee  s  (1;  2  =  Eee+  (3-2) 

where  E  is  the  sign  of  expectation  value,  and  sign  +  refers  to  the  matrix  or 
vector  transposition. 

At  present,  the  numerous  methods  to  find  a  stable  solution  of  the  Fredholm 
equation  of  the  first  kind  (3.1)  have  been  developed.  One  of  the  commonly  used 
inversion  method  is  the  statistic  regularization  techniques  (method  of,  optimal 
estimation)  proposed  independently  by  several  authors  (Strand  and  Westwater,  1968; 

Turchin  et  al.,  1968:  Turchin  and  Nozik,  1969;  Rodgers,  1970,  1971,  1976,  1990). 
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To  solve  ill-posed  problem  (3.1)  requires  regularizing  the  equation  on  the  basis 
of  additional  information  about  decision. 

In  the  linear  statistical  method  of  the  regularization,  the  decision  is  assumed  to 
belong  to  an  a  priori  statistics,  defining  the  average  vector  of  physical  state 
parameters  x,  and  its  covariance  matrix  K  describing  a  natural  variability  of  the 
sought  parameter  profiles: 

x  =  Ex;  K  =  E(x  -  x)(x  -  x)+  (3.3) 

The  covariance  matrix  K  defines  the  valuation  ||x  -  xl^  characterizing  a 
departure  of  vector  x  from  the  statistics  (3.3): 

II x  -  *|i-i  =  (x  -  x,  K'x(x  -  x))  (3.4) 

lx 

Then  the  estimate  x  of  regularized  problem  (3.1)  can  be  obtained  by  joint 
minimizing  the  least-square  fit  to  the  measurements  y  and  the  least-square  fit  to 
the  statistical  mean  x: 

II  y  -  Ax  1^,-1  +  ||x  -  xf^  =>  min  (3.5) 

or,  what  is  the  same, 

(y  -  AxdtT^y  -  Ax))  +  (x,K'xx)  =>  min  (3.6) 

(For  simplicity  in  (3.6)  and  further  we  have  neglected  by  the  shift  transformation 
of  X  and  Y  spaces,  i.e.  substituted,  for  example,  x  -  x  by  x). 

Minimizing  (3.6)  with  respect  to  x  yields: 

Gx  +  K  1x  =  A+2  *y  (3-7) 

where  G  =  A+2  *A  is  the  Fisher  matrix. 

From  (3.7)  one  can  receive  the  following  expressions  for  the  retrieval  estimate  x 
and  its  covariance  error  matrix  D: 

x  =  x  +  (G  +  K’VVS'V  (3.8) 

D  =  (G  +  K'1)'1 

To  estimate  the  accuracy  of  the  linear  solution  and  the  information  content  of  the 
measurements,  V.P. Kozlov  (1966)  have  brought  into  existence  the  matrix  C,  named  as 
"information  matrix": 

C  =  KG  (3.9) 
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If  Uj  and  X;  are  the  eigenvectors  and  eigenvalues  of  the  matrix  C  (Cu=Xu),  then 

m 

the  coefficients  b;  in  expansion  of  vector  x  =  £  bjU;  are  incorrelated  with  each 

i=i 

other  and  independently  estimated  values. 

There  exist  a  variety  of  approaches  to  estimate  the  information  content  of  the 
measurements  (Westwater  and  Strand,  1970;  Houghton  et  al.,  1986).  We  have  used  two 
definitions  of  information  content  proposed  by  Shannon  (Shannon  and  Weaver,  1949) 
and  Kozlov  (Kozlov,  1966). 

The  information  content,  according  to  Shannon,  may  be  regarded  roughly  as  the 
number  of  bits  needed  to  represent  the  number  of  distinct  measurements  that  could 
have  been  made.  This  is  computed  as  (Houghton  et  al.,  1986): 
m 

I  =  0.5  [log2(l+Xi)  .  (3.10) 

1 

The  "information  volume",  according  to  Kozlov,  is  defined  as  (Kozlov,  1966): 

V  = 

Evidently,  when  calculating  the  information  volume,  we  use  only  eigenvalues  which 
are  greater  than  1  and  correspond  to  the  information  components  distinguished  at 
given  level  of  the  measurement  error.  The  point  is  that  each  atmospheric  remote 
sensing  measurement  is  practically  unique,  and  statistical  selection  of  the  signal 
components  which  are  smaller,  than  measurement  error,  is  impossible.  And 
consequently,  the  information  volume,  according  to  Kozlov,  may  be  regarded  as  the 
number  of  atmospheric  states  distinguishable  by  used  measurement  device  and  may 
reach  to  very  large  magnitudes. 

Another  usefull  measure  of  information  content  (named  as  "Kozlov  factor")  is  the 
number  of  eigenvalues  exceeding  in  magnitude  1.  This  factor  is  the  quantity  of 
independent  information  components  that  could  have  been  derived  from  the 
measurements  performed  with  given  error. 

It  is  evident,  that  in  the  case  of  uncorrelated  measurement  error  (i.e.  if  2  = 

a2 1,  where  I  is  the  unit  matrix),  the  above  eigenvalues  X;  and  eigenvectors  u;  are 

connected  with  eigenvalues  p;  and  vectors  <j>[  of  the  covariance  matrix  of  measured 

quantity  K  =  E(yy+)  by  the  simple  relations: 

X;  =  [i; /a2  <(>;  =  Yl\u-  (3.12) 

(In  discussed  case  K  =  E(yy+)  is  K  =  E(TbtJ).) 

tbtb 
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The  above  expressions  for  the  information  content  estimations  have  been  defined 
for  the  case  of  linear  operator  A.  As  a  rule,  in  the  atmospheric  applications,  the 
problem  is  described  by  linear  approach  only  with  some  error.  Nevertheless,  in 
view  of  Eq.(3.12),  we  can  use  K  matrix  calculated  with  regard  for  nonlinear 
relations  between  outgoing  radiance  and  atmospheric  parameters  to  estimate  the 
information  content  /  and  information  volume  V. 

3.2.The  basis  for  the  calculations  of  outgoing  microwave  radiance 
information  content 

\ 

The  radiative  transfer  model  described  in  the  Chapter  2  was  used  for  calculations 

of  covariance  matrices  of  outgoing  radiance  in  the  terms  of  brightness  temperature 

(K  )  in  the  1-200  GHz  spectral  range. 
tbtb 

To  calculate  K  ,  it  is  necessary  to  have  some  simulated  set  of  the  atmospheric 
tbtb 

states  correctly  describing  the  real  atmospheric  and  surface  variabilities. 

In  particular,  an  empirical  ensemble  corresponding  to  a  test  area  can  be  used  as 
that  set.  One  can  also  carry  out  a  statistical  simulation  of  the  set  using  the 
values  of  average  profiles  and  covariance  matrices,  generalizing  the  long-term 
seasonal  observations  in  the  test  region. 

Another  method  for  forming  the  standard  set  has  been  proposed  and  successfully 

used  in  our  work  (Kozlov  et  al.,  1983).  This  approach  uses  some  ideas  of  classical 

theory  of  planning  the  experiment,  in  particular,  the  well-known  Box  theorem  (Box, 

1952).  In  the  present  study,  to  calculate  K  we  have  utilized  a  priori  set, 

tbtb 

formed  by  Kozlov  et  al.  (1983)  from  "global"  covariance  matrices  of  temperature 
and  humidity,  constructed  by  Dr.  Arking  (Goddard  Center  of  Space  Investigation). 
The  cross-correlations  of  these  parameters  have  assumed  to  be  equal  to  0,  that 
extends  a  statistical  distribution. 

As  a  result,  we  received  the  standard  a  priori  set,  consisting  of  40  atmospheric 

models  (temperature,  humidity  profiles)  for  quasi-statistical  method.  This  set 

describes  a  priori  statistics  of  the  atmospheric  state  variations  in  a  generally 

correct  way  (relatively  to  the  first  two  moments  of  the  distribution). 

To  estimate  a  range  of  applied  atmospheric  state  ensemble,  it  is  sufficient  to 

mention,  that  the  mean  square  variations  of  atmospheric  temperature  and  humidity 

at  different  levels  reach  to  16.6K  and  80%,  respectively,  and  corresponding  value 

for  surface  temperature  variations  is  16.9K.  So,  applied  set  is  "very  wide". 

For  the  information  content  analysis,  K  calculations  were  carried  out  for 
.  tbtb 
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monochromatic  case  and  1500  frequencies  (so-called  "the  complete  radiance 
experiment").  Their  locations  and  amount  in  different  intervals  are  given  in  Table 
3.1. 


Table  3.1 

Frequency  assignment  for  outgoing  radiance  calculations 


Spectral  range 
GHz 

Absorber 

Number  of  the 
points 

40-85 

02,h20 

575 

110  -130 

h20 

175 

15  -  30 

h2o 

125 

175  -195 

h2o 

125 

'  The  additional  500  measurement  channels  are  distributed  uniformly  in  all  spectral 
range  of  1-200  GHz.  ■ 

As  it  follows  from  cited  table,  chosen  frequencies  cover  oxygen  absorption  range, 
transparency  windows  and  water  vapor  absorption  lines. 
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3.3. Estimations  of  the  number  of  independent  information  components 

In  Table  3.2  we  give  the  first  30  eigenvalues  p:  of  K  matrix  calculated  with 

tbtb 

regard  for  atmospheric  temperature  and  humidity  variations  and  surface  emissivity 
equal  to  0.95. 


The  first  eigenvalues  pj  (K‘)  of  K 


Tbtb 


Table  3.2 

matrix  for  the  1-200  GHz  range. 


(1500  channels;  40  different  temperature  and  humidity 
profiles  and  surface  temperatures;  land,  e  =  0.95) 


1 

256467.69059 

2 

9458.45364 

3 

8967.11645 

4 

647.28262 

5 

546.44578 

6 

467.17321 

7 

264.37950 

8 

200.09120 

9 

91.15060 

10 

38.50167 

11 

13.75746 

12 

7.54997 

13 

3.10991 

14 

2.12617 

15 

1.61347 

16 

0.81480 

17 

0.47060 

18 

0.31865 

19 

0.24341 

20 

0.17757 

21 

0.12955 

22 

0.09117 

23 

0.07151 

24 

0.03988 

25 

0.03029 

26 

0.01810 

27 

0.01530 

28 

0.01178 

29 

0.01013 

30 

0.00644 

It  follows  from  these  data,  that  15  p;  (for  aT  =1.0K),  18  Pj  (for  aT =0.5K)  and  29 

AB 

p;  (for  oTb=0.1K)  comply  with  V.P. Kozlov  criterion: 

h  >  °tb  • 

where  OjQ  is  the  random  noise  variance  of  outgoing  radiance  measurement  in  the 


42 


terms  of  TB. 

Similar  Kt  t  and  p,;  calculations  for  linear  spectral  model  of  surface  emissivity 

of  the  sea  (e(f)=a  +  bf)  result  in  the  following  number  of  independent  components: 
18  (for  oTb=1.0K),  22  (for  ctTb=0.5K),  and  31  (for  0^=0. IK). 

So  the  1500  monochromatic  spectral  measurements  simulating  the  complete  radiance 
experiment  contain  a  great  deal  of  information  about  vertical  temperature  and 
humidity  profiles  of  the  atmosphere  and  underlying  surface  temperature  Ts.  It  is 
worth  to  emphasize,  that  these  calculations  did  not  take  into  account  the  possible 
variations  of  surface  emissivity. 

The  conducted  calculations  do  not  allow  to  separate  the  "information  components" 

of  matrix  K  associated  with  the  variations  of  different  atmospheric  parameters 
tbtb 

(T(z)  and  q(z)). 

In  this  connection,  we  carried  out  the  additional  K  computations  when  varying 

tbtb 

T(z)  (at  q(z)=q(z))  and  q(z)  (at  T(z)=T(z))  in  turn.  The  corresponding  eigenvalues 
(at  surface  emissivity  equal  to  0.95)  are  given  in  Table  3.3  and  3.4: 


i 


/ 
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Table  3.3 


The  first  eigenvalues  p,;  (K‘)  of  K  matrix  for  the  1-200  GHz  range. 


tbtb 


(1500  channels;  40  different  atmospheric  temperature  profiles  and  surface 
temperatures;  the  mean  atmospheric  humidity  profile;  land,  e  =  0.95) 


1 

254144.99630 

2 

8975.13144 

3 

2102.01714 

4 

547.25454 

5 

210.81622 

6 

127.15609 

7 

38.61295 

8 

13.68913 

9 

6.68457 

10 

2.25899 

11 

1.39076 

12 

0.68921 

13 

0.56537 

14 

0.32975 

15 

0.19811 

16 

0.09753 

17 

0.03694 

18 

0.02981 

19 

0.01202 

20 

0.00673 

7 

The  first  eigenvalues  p;  (K")  of  K 


tbtb 


Table  3.4 

matrix  for  the  1-200  GHz  range. 


(1500  channels;  40  different  atmospheric  humidity  profiles; 
fixed  surface  temperature  and  the  mean  atmospheric 
temperature  profile;  land,  e  =  0.95) 


7639.21629 

820.94697 

281.12041 

68.05590 

5.73913 

0.33495 

7 

0.10836 

8 

0.04536 

HI 

0.01694 

m 

0.01133 

H 

0.00672 
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The  analysis  of  table  data  has  showed  the  following: 

1. The  numbers  of  independent  information  components,  inherent  in  outgoing  radiance 
measurements  over  land  and  associated  with  only  atmospheric  and  surface 
temperature  variations,  are  11,  14,  and  19  at  levels  of  measurement  errors  (ctTb) 
are  1.0K,  0.5K  and  0.1K,  respectively. 

2. The  numbers  of  independent  information  components,  inherent  in  outgoing  radiance 
measurements  over  land  and  associated  with  only  atmospheric  humidity  profile 
variations,  are  5,  6,  and  10  at  equal  to  1.0K;  0.5K  and  0.1K,  respectively. 

3. Replacing  the  land  emissivity  model  by  calm  sea  one  only  slightly  increases  (by 
1-2)  the  number  of  distinguished  independent  components  in  the  outgoing  radiance 
spectrum. 

4. The  sum  of  independent  component  numbers  of  the  covariance  matrices  calculated 

with  regard  for  T(z)  or  q(z)  variations  individually  is  equal  approximately  to  the 

number  of  those  for  K  covariance  matrix  considering  these  variations  in 

tbtb 

combination.  Apparently,  the  last  is  a  result  of  neglecting  the-  interconnection 
between  atmospheric  temperature  and  humidity  when  composing  the  a  priori  set  of 
atmospheric  state  ensembles. 

The  principal  conclusion  from  these  first  stage  studies  can  be  formulated  in  the 
following  way.  The  complete  radiance  experiment  contains  a  considerable  amount  of 
information  about  T(z)  and  q(z)  profiles  in  terms  of  the  principal  components  of 
the  covariance  matrix  of  outgoing  brightness  temperatures.  Their  number  exceeds 
noticeably  the  number  of  SSM/T  and  SSM/T-2  measurement  channels  (12  channels).The 
number  of  independent  components,  which  could  be  potential-derived  from  outgoing 
radiance  spectra,  has  to  make  possible,  in  principle,  to  retrieve  temperature  and 
humidity  profiles  with  high  accuracy  provided  a  low  level  of  measurement  noise 
(ctTb=0.1K).  Table  3.5  summarizes  the  data  about  the  number  of  the  independent 
components  for  different  calculation  versions  and  error  levels. 


Table  3.5 

Number  of  independent  components  at  different  calculation  variants 


No. 

Calculation  variant 

°TB-  k 

0.1 

0.5 

1.0 

1. 

T(z),  Ts  and  q(z) 
variations,  e=0.95 

29 

18 

15 

2. 

T(z),  Ts  and  q(z) 
variations,  e(f)=a+bf 

31 

22 

18 

3. 

T(z)  and  T, 
variations,  e=0.95 

19 

14 

11 

4. 

T(z)  and  Ts 
variations,  e(f)=a+bf 

19 

15 

12 

5. 

q(z)  variations, 
e=0.95 

10 

6 

5 

6. 

q(z)  variations, 
e(f)=a+bf 

11 

7 

6 

It  is  interesting  to  compare  the  information  content  of  the  complete  radiance 

experiment  with  that  of  different  specific  satellite  devices.  For  the  correctness 
of  comparison,  the  measurements  by  different  devices  were  simulated  without 

considering  their  finite  spectral  resolutions  and  for  the  monochromatic  case.  In 
our  calculations,  we  used  the  means  of  monochromatic  values  in  the  centers  of 

corresponding  intervals  for  a  measurement  channel  consisting  of  two  or  more 
spectral  sub-intervals.  The  special  study  of  this  question  has  shown,  that  using 
the  "monochromatic  approach"  for  examined  satellite  devices  do  not  bring  the 

noticeable  errors  in  calculations  of  the  informative  content. 

In  Table  3.6  the  first  eigenvalues  of  K  matrix,  calculated  using  the  above 

tbtb 

a  priori  atmospheric  set,  are  given  for  a  number  of  satellite  devices. 
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Table  3.6 


a)  Land 


b)  Sea 


2 

The  first  eigenvalues  |x;  (K )  of  K  matrix  for 

tbtb 

a  number  of  satellite  devices. 


Devices 

i 

SSM/T 

SSM/T-2 

SSM/T 

+ 

SSM/T-2 

1 

347.58 

949.03 

1297.23 

2 

175.79 

127.87 

179.11 

3 

8.40 

11.30 

149.27 

4 

4.25 

1.03 

10.39 

5 

0.58 

0.45 

8.78 

6 

0.15 

5.87 

7 

0.0074 

1.19 

8 

0.94 

9 

0.52 

10 

0.13 

11 

0.019 

12 

0.0034 

Devices 

i 

SSM/T 

SSM/T-2 

SSM/T 

+ 

SSM/T-2 

1 

215.60 

743.98 

910.38 

2 

174.70 

560.72 

612.38 

3 

18.24 

26.17 

178.17 

4 

7.96 

4.30 

33.66 

5 

1.15 

0.68 

8.21 

6 

0.24 

5.26 

7 

0.025 

3.37 

8 

1.22 

9 

0.61 

10 

0.38 

11 

0.099 

12 

0.0064 

An  analysis  of  data  of  the  table  3.6  let  us  make  the  following  conclusions: 
l.For  radiometer  SSM/T,  the  numbers  of  distinguished  (according  to  Kozlov 
criterion)  independent  components  in  the  covariance  outgoing  radiance  matrix  make 
up  4,  5  and  6  (for  land)  and  5,  5  and  7  for  sea  at  RMS  random  errors  cr^of 
brightness  temperature  measurement  equal  to  1.0,  0.5  and  0.1K,  respectively. 

Accounting  for  the  vicinity  of  seventh  eigenvalue  for  land  to  10-2  (0.0074),  it 
can  be  concluded,  that  all  SSM/T  channels  at  measurement  error  0.1K  are 
"informative"  ones.  At  higher  error  levels  (1.0  and  0.5K),  radiometer  SSM/T  has 
the  "redundant"  measurement  channels.  Having'  in  mind  the  values  of  informative 
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components  given  in  table  3.5  (15,  18  and  29  for  land),  we  can  conclude  that  this 
fact  testifies  to  non-optimal  spectral  positions  of  SSM/T  channels. 

2.For  radiometer  SSM/T-2,  the  numbers  of  informative  components  are  4,  5  and  5 
(over  sea  and  land)  at  error  levels  equal  to  1.0,  0.5  and  0.1K.  So,  even  at 
ctTb=0.5K,  all  SSM/T-2  channels  are  informative. 

For  a  case  of  joint  use  of  SSM/T  and  SSM/T-2  radiometers,  the  number  of 
independent  components  is  given  in  table  3.7. 

Table  3.7 

The  number  of  independent  components  of  outgoing  radiance 
covariance  matrix  at  joint  use  of  SSM/T  and  SSM/T-2 

radiometers 


Measurement  error,  K 

0.1 

0.5 

1.0 

Land 

11 

9 

7 

Sea 

11 

10 

8 

As  it  follows  from  table  3.7,  for  combined  remote  sensing  system  of  DMSP  satellite 
(SSM/T  and  SSM/T-2),  all  the  measurement  channels  are  practically  informative  at 
0^=0. IK.  At  higher  error  levels  (0.5  and  1.0K),  the  channels  of  this  remote 
sensing  system  are  not  all  informative.  This  fact,  with  accounting  for  data  of 
table  3.5,  testifies  indirectly  to  non-optimal  choice  of  the  channels  for 
temperature-humidity  sounding  by  DMSP  satellite.  It  is  worth  note,  that  these 
conclusions  are  concerned  to  examined  "global  ensemble"  of  atmospheric  states  and 
the  cloudless  atmosphere.  Certainly,  the  inclusion  of  the  variations  of  cloud  and 
surface  characteristics  can  influence  this  result. 

3.4.Information  characteristics  of  the  measurements  by  different  devices 

Using  calculated  eigenvalues  of  outgoing  radiance  covariance  matrix  and  the  Eqs. 
(3.10)-(3.12),  we  have  determined  the  information  content  I  (by  Shannon)  and 
information  volume  V  (by  Kozlov)  of  different  outgoing  radiance  measurements.  The 
examples  of  these  calculations  at  measurement  RMS  error  of  0.1K  are  given  in  table 
3.8. 

(In  the  tables  3.8,  3.9  letter  E  denotes  the  following:  for  example,  0.346E28  = 
0.346- 10+28). 
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Table  3.8 

Information  content  I  (according  to  Shannon)  and  information 
volume  V  (according  to  Kozlov)  for  different  spectral  designs 
of  the  measurements  at  aT  =0.5K 

■*B 


Land 

Sea 

scheme  of 
the  measurements 

I  V 

I  V 

1500  channels, 

all  parameters  are  variable 

74.2  0.249E+22 

86.0  0.103E+26 

1500  channels, 
temperature  is  variable 

54.2  0.519E+16 

54.8  0.819E+16 

1500  channels, 
humidity  is  variable 

25.8  0.307E+08 

33.8  0.866E+10 

SSM/T 

all  parameters  are  variable 

15.9  0.375E+05 

17.0  0.801E+05 

SSM/T-2 

all  parameters  are  variable 

15.1  0.256E+05 

17.7  0.181E+06 

SSM/T  +  SSM/T-2 

all  parameters  are  variable 

26.7  0.532E+08 

29.6  0.383E+09 

From  these  data  it  can  be  drawn  the  following  conclusions: 

l.Information  content  of  the  "complete  radiance  experiment"  (1500  measurements  in 
the  1-200  GHz  spectral  range)  is  essentially  higher  than  those  of  SSM/T  and  SSM/T- 
2  radiometer  measurements  (both  using  those  individually  and  in  combination). 
2.Information  content  of  outgoing  radiance  measurements  over  water  surface  with 
respect  to  vertical  atmospheric  and  humidity  profiles  is  higher  than  that  over 
land. 

The  information  data  given  in  table  3.8  does  not  account  for  essentially  different 
number  of  outgoing  radiance  measurements  in  the  complete  radiative  experiment 
(N1=1500)  and  for  SSM/T  and  SSM/T-2  radiometers  (N2).  But  it  is  evident,  that 
increasing  a  number  of  the  measurements  (even  in  the  same  spectral  channels) 
decreases  practically  a  level  of  random  measurement  errors  (in  relation  of  Nj/N2' 
for  the  random  error  variance).  In  this  connection,  we  give  table  3.9  representing 
the  measurement  information  content  according  to  Shannon  and  Kozlov  with 
accounting  for  different  number  of  the  measurements  in  various  experiment  designs. 
In  the  Table  3.9,  cXb  for  different  devices  are  calculated  from  oTb(N2)  = 


On 


a(Nx) 


N, 


and  correspond  to  aTB(Nj)=0.4K. 
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Table  3.9 

Information  content  /  and  information  volume  V 
for  different  spectral  designs  of  measurements. 

(A  priori  ctTb  for  1500  channels  -  0.4K,  for  other  designs 

°tb  *s  given  *n  the  first  column  of  the  table) 


Land 

Sea 

scheme  of 
the  measurements 

I  V 

I  V 

1500  channels, 

all  parameters  are  variable 

81.4  0.343E+24 

94.2  0.301E+28 

1500  channels, 
temperature  is  variable 

59.4  0.213E+18 

60.1  0.378E+18 

1500  channels, 
humidity  is  variable 

28.1  0.142E+09 

36.5  0.517E+11 

SSM/T 

all  parameters  are  variable 
aT  =  0.026  K 

42.0  0.427E+13 

44.3  0.210E+14 

SSM/T-2 

all  parameters  are  variable 
aT  =  0.024  K 

37.0  0.143E+12 

39.9  0.101E+13 

SSM/T  +  SSM/T-2 

all  parameters  are  variable 

ctTb  =  0.035  K 

66.7  0.102E+21 

72.6  0.317E+22 

Table  3.9  demonstrates,  that  information  content  is  noticeably  higher  in  the 
"complete  radiance  experiment"  than  in  the  DMSP  temperature-humidity  sounding 
system.  For  example,  information  contents  by  Shannon  are  81.4  (over  land),  94.2 
(over  sea)  and  66.7  (over  land),  72.6  (over  sea)  for  the  "complete  radiance 
experiment"  and  for  simultaneous  SSM/T+SSM/T-2  measurements,  respectively. 

So,  the  principal  result  of  performed  studies  consists  in  revealing  the  fact,  that 
DMSP  temperature-humidity  sensing  system  does  not  make  possible  to  receive  all 
available  information  about  atmospheric  '  temperature  and  humidity  profiles 
contained  in  outgoing  microwave  radiation  spectra. 

In  conclusion  it  is  necessary  to  note,  that  we  have  carried  out  the  similar 
information  content  analysis  for  20-channel  AMSU  radiometer  to  be  launched  to 
space  flight  in  1996.  The  calculations  have  shown  the  noticeable  preference  the 
information  characteristics  of  AMSU  in  comparison  with  those  of  SSM/T+SSM/T-2 
system.  But  the  information  content  of  AMSU  measurements  is  appreciably  smaller 
than  that  of  the  "complete  radiance  experiment". 


50 


4. ANALYSIS  OF  POTENTIAL  ACCURACY  OF  TEMPERATURE  AND 
HUMIDITY  MICROWAVE  REMOTE  SOUNDING 

4.I.Mathematical  basis  of  estimating  the  potential  accuracy  of 
remote  sensing 

In  the  case'  of  linear  operator  A  (see  Part  3.1),  the  solution  of  the  inverse 
problem  may  be  represented  according  to  (3.8)  as: 

x  =  x  +  (G  +  K'YA^y 


Then  we  can  write  the  following  expression  for  the  error  of  inverse  solution: 

Sx  =((G  +  K'1)'^  -  I)x  +  (G  +  K‘1)‘1A+2"1e  (4.1) 

here  G  =  A+2_1A  is  the  Fisher  information  matrix,  K  is  the  covariance  matrix  of 
the  a  priori  variability  of  atmospheric  parameter  vector  x,  I  is  the  unit  matrix, 
2  is  the  noise  covariance  matrix. 

It  is  easily  seen,  that  the  first  term  in  Eq.4.1  does  not  depend  on  concrete 
realization  of  the  random  measurement  error  e,  and  the  second  -  on  the  sought 

vector  of  atmospheric  parameters  x.  Therefore,  it  is  the  practice  to  call  the 

first  term  as  systematic,  or  null-space  error  component  (Rodgers,  1990),  and  the 

second  term  -  as  random  component.  As  those  are  independent  random  values,  the 
covariance  error  estimate  matrix  D  (see  formula  (3.9))  may  be  represented  as  the 

sum  of  their  covariance  matrices: 


D  -  Dr  +  Ds 

(4.2) 

DR  =  (G  +  K“1)'1G(G  +  K'1)'1 

(4.3) 

Ds  =  ((G  +  K'1)'^  -  I)K((G  +  KVG  -  I) 

(4.4) 

where  DR  and  Ds  are  the  covariance  matrices  of  the  random  and  systematic  retrieval 
error  components,  respectively. 

The  expressions  (3.8),  (4.1)  have  been  formulated  for  one-vector  retrieval  problem 
(i=l)  (Strand  and  Westwater,  1968;  Turchin  and  Nozik,  1969;  Turchin  et  al.,  1970; 
Rodgers,  1970,  1971,  1976).  But  in  reality  the  retrieval  problems  are  many- 
parametric  ones,  and  separate  studies  of  different  retrieval  problems  present,  as 
a  rule,  only  a  methodical  interest,  for  example,  to  examine  the  maximum 
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possibilities  of  the  retrieval  methods.  But  even  in  these  separate  retrievals  it 
is  necessary  to  take  into  account  not  only  measurement  errors  (noise  errors)  but 
"effective"  noise  due  to  the  variations  of  other  parameters,  though  the 
retrievals  of  those  are  not  planned. 

We  will  illustrate  this  fact  by  formulation  of  the  unified  atmospheric  temperature 
(T)  and  humidity  (q)  profile  retrievals  (Pokrovsky  and  Timofeyev,  1972),  because 
this  problem  relates  just  to  microwave  sensing  of  the  cloudless  atmosphere  by 
spectrometers  SSM/T-1,  SSM/T-2.  In  this  case,  assuming  a  known  emissivity  of 
underlying  surface  we  can  formulate  the  inverse  problem  in  the  following  way 
instead  of  Eq.(3.1): 

by  =  At  6T  +  Aq  6q  +  e  (4.5) 


Then  from  (3.8)  one  can  obtain  for  optimal  statistical  estimates  of  temperature, 
6T,  and  humidity,  bq,  variances  : 

~  +  +  +  +  + 

6T  =  (Kyy  Ay  +  Krpq  A^)  (Ay  ^TT  Ay  +  Ay  ^Tq  Aq  +  A^  K^y  Ay 

+  Aq  Kqq  Aq  +  S)"1  by  (4.6) 

„  +  +  +  +  + 

bq  =  (Kqq  Aq  +  KqT  At)  (At  Kyp  AT  +  AT  KTq  Aq  +  Aq  KqT  AT  + 

+  A,  K„  +  S)-‘  6y  (4.7) 


where  Kyp,  Kqq  and  are  covariance  and  cross -covariance  matrices  of  the 

sought  parameters. 

Analogously  from  (3.9)  we  can  deduce  for  error  matrices  characterizing  the 
accuracies  of  the  retrieval  methods: 


Dn  -  ] 


-l 


-l 


Krp  +(At  +  K-pp  KTq  Aq)  [Aq(Kqq  -  KqT  Kyp  KTq)Aq  + 


+  Z]'1  (AT  +  Aq  KqT  K^) ' 


(4.8) 


K  +  Kl  KqT  A^)  [A^Kyp  -  KTq  Kl  KqT)Ay  + 


+  Z]'1  (A,  +  AT  KTq 


(4.9) 
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Supposing  Aq  =  0  and  KTq  =  0,  we  can  write  the  formulas  (4.6),  (4.8)  for 
temperature  profile  retrieval  (one-dimensional,  idealized  problem): 


~  + 

6T  =  K yy  At  (At  Kjj 


At 


+  J)'1  6y 


(4.10) 


D'it  - 


yl 

'  Kjp  +  AT2  At  ’ 


(4.11) 


Comparing,  for  example,  (4.6)  and  (4.10)  one  can  realize  that  idealized  approach 
suffers  from: 

-  distortion  of  the  operator  of  the  inverse  problem; 

-  lack  of  the  term  Aq(Kqq  -  KqT  K-^  KTq)Aq  presenting  by  itself  an  additional 
"effective"  noise  in  the  radiance  measurements. 

In  our  study  we  assume,  that  temperature  and  humidity  profiles  are  independent, 
that  exceeds  the  a  priori  statistical  ensemble.  This  means  that  KqT  =  KTq  =  0. 
Then  Eqs.(4.8)  and  (4.9)  turn  into: 


Dtt  - 


K-rr  +  AT(AqKqqAq  +  2)  AT 

k. 

*  1 
Kqq  +  Aq(ATKTrAT  +  2)  Aq 


(4.8') 

(4-9') 


From  these  equations  it  is  evident,  that  the  uncertainty  of  one  atmospheric 
parameter  influence  on  retrieval  error  of  other  parameter.  In  Eq.(4.8),  for 

example,  the  term  A^Kqq-KqjKypKj^A^,  added  to  covariance  measurement  matrix, 
represents  by  itself  an  "effective"  noise  contaminating  temperature  profile 
measurements  because  of  the  humidity  profile  uncertainties.  In  the  case  of 
neglecting  the  cross-correlations  between  temperature  and  humidity  (4.8 ),  this 

+ 

effective  noise  is  AaK_aAn. 

The  above  studies  of  device  measurement  content  (Chapter  3)  and  retrieval  error 
have  been  based  on  the  assumptions  of  linear  relations  between  radiance  and  the 
sought  atmospheric  parameters.  But  the  basic  retrieval  problems  of  the  remote 
sounding  are  nonlinear.  For  example,  in  the  1-200  GHz  spectral  range,  the 
dependence  between  outgoing  brightness  temperature  and  atmospheric  humidity 
profile  is  essentially  nonlinear.  Formally  the  nonlinear  problem  can  be  written  as 
a  generalization  of  the  linear  problem  (3.1): 
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y  =  B(x)  +  e 


(4.12) 


where  B  is  nonlinear  operator  and 

x  =  Ex;  e  =  Ee  =  0;  K  =  E(x  -  x)(x  -  x)+;  2  =  Eee+ 

And  the  best  estimate  (in  the  probability  sense)  x  corresponding  to  measurement 
value  y  minimizes  the  sum  of  the  following  deviations: 

|B<x)  -  yl’.,  +  | x  -  =»  min  (4.13) 

The  traditional  way  to  solve  a  nonlinear  problem  is  to  reduce  the  task  to  the 
sequence  of  linear  tasks  by  some  form  of  the  iteration.  The  iteration  method  is 
simply  a  matter  of  expanding  the  operator  B  as  a  Taylor  series  about  a  guessed 
value  xk: 

B(x)  =  B(xk)  +  A(x  -  x,J  +  o(x  -  x^ 

here  A  is  variational  derivative  matrix  of  operator  B  at  the  point  xk,  xk  is  the 
solution  estimate  derived  on  the  previous  iteration  step. 

Then  Eq.(4.13)  reduces  to: 

||B(xk)+A(Xk+1-Xk)-y||‘,-1  +  |xk+1-x  |  *  =>  min  (4.14) 

2  K  \ 

By  minimizing  this  equation  one  can  derives  some  equivalent  and  useful  in  practice 
expressions  for  xk+1: 


xM=(A+r'A  +  K-y‘(A+r'(y  -  B(Xk)  +  Axt)  +  K*1*)  (4.15-1) 

xlu=;  +  (AV'A  +  K-y‘A+r‘(y  -  B(Xk)  -  A(i-Xk))  (4.15-2) 

xttI=x1+(A+r1A  +  K-V'(AV(y  -  B(xk))+K"1(x-xlt))  (4.15-3) 

xklI=i  +  KA+(AKA+  +  Sj'Vy  -  B(Xk)  -  Alx-x,))  (4.15-4) 


It  is  easy  to  prove  that  the  limit  of  decision  sequence  (if  the  process  converges) 
is  the  solution  of  nonlinear  problem  (4.13). 

Although  Eqs.  (4.15-l)-(4.15-4)  are  equivalent,  the  proper  choice  of  applied 
formula  can  essentially  simplify  a  development  of  relevant  software.lt  should  be 
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mentioned  that  similar  iterative  procedures  were  early  developed  by  other  means 
(for  example,  Marks  and  Rodgers,  1993;  Polyakov  and  Rozanov,  1989). 

In  a  number  of  cases,  the  above  algorithm  does  not  give  a  convergent  sequence  of 
{xk}  due  to  essential  nonlinearity  of  the  operator  B.  To  resolve  the  problem, 
authors  of  the  paper  (Marks  and  Rodgers,  1993)  have  proposed  to  use  the  known 
Marquard  algorithm  of  the  functional  minimizing,  that  gives  the  following 
iterative  scheme: 

xk+1=xk+(A+2"1A+K'1+cI)'1(A+2'1(y  -  B(xk))+K‘1(x-xk))  (4.16) 

where  c  is  an  adjusting  coefficient,  I  is  the  unit  matrix. 

We  have  proposed  (Polyakov,  1995)  the  approach  somewhat  different  from  the  above 
(and  stated  further)  to  solve  this  problem.  It  gives  the  expression  (4.18), 
similar  to  (4.16),  but  having  a  simple  physical  sense  and  some  advantages  in 
using.  This  approach  is  the  following. 

In  the  case  of  iterative  process  divergence,  it  is  the  most  reasonable  to  restrict 
the  solution  variations  at  each  iterative  step  by  means  of  adding  the  limiting 
term  to  minimize  the  expression  (4.14),  i.e.: 

II  B(xk+i)-y  II |,-i+ 1  Xk+1-X  1  ^.-1+  II  xk+1-xk  ||  =>  min  (4.17) 

Here  L  is  some  positive  definite  matrix.  Needless  to  say,  Eq.(4.17)  is  miningful 
only  in  iterative  process. 

If  the  iterative  process  converges,  the  last  term  in  Eq.(4.17)  tends  to  •  zero, 
whatever  matrix  we  use.  This  matrix  may  even  change  from  iteration  step  to  step. 
In  any  case,  the  solution  of  (4.17)  coincides  with  that  of  the  problem  (4.14). 

Then  the  formula  of  the  solution  refinement  similar  to  (4.15-3)  is: 

xk+i=xk+(A+2'1A+K’1+L~1)"1(A+2"1(y  -  B(xk))+K'1(x-xk))  (4.18) 

By  analogy  with  matrix  K  characterizing  the  estimate  deviation  from  the  mean, 
matrix  L  is  the  covariance  matrix  of  the  ^estimate  variation  at  each  iterative 
step.  Matrix  L  must  be  constructed  with  regard  for  natural  atmospheric  parameter 
variability  and  the  degree  of  nonlinearity  in  the  radiance-parameter  dependencies. 
For  example,  in  the  examined  task,  the  dependence  between  radiance  and  atmospheric 
temperature  profile  is  close  to  linear,  but  that  between  radiance  and  atmospheric 
humidity  profile  is  essentially  nonlinear.  So,  matrix  L  may  be  a  diagonal  one  with 
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large  and  small  components,  corresponding  to  temperature  and  humidity, 
respectively.  In  each  specific  case,  it  is  appropriate  to  construct  this  matrix  on 
the  basis  of  special  feature  of  the  task  and  to  fit  its  optimal  components  by 
numerical  simulation. 

4.2.Numerical  estimation  of  the  potential  errors  of  the  temperature 
and  humidity  retrievals 

The  above  study  of  microwave  measurement  information  content,  presented  in  Chapter 
3,  has  demonstrated  an  availability  of  additional  information  about  T(p)  and  q(p) 
in  the  radiance  compared  to  DMSP  remote  sensing  measurements  (SSM/T  and  SSM/T-2 
radiometers). 

To  estimate  quantitatively  an  accuracy  of  atmospheric  temperature  and  humidity 
sounding,  we  have  carried  out  the  error  matrix  calculations  using  the  basic 
statistical  approach  (Pokrovsky  and  Timofeyev,  1972).  To  this  end,  two  designs  of 
outgoing  microwave  radiance  measurements  were  examined: 

a) ."Complete  radiance  experiment"  -  1500  measurement  channels  in  the  whole 
spectral  range  of  1-200  GHz  (see  Table  3.1); 

b) . SSM/T  and  SSM/T-2  measurements  (12  measurement  channels). 

The  large  volume  of  error  matrix  calculations  has  made  possible  to  compare  two 

experiment  designs,  to  study  the  effects  of  random  measurement  errors,  vertical 

resolution,  and  a  priori  statistics.  Basic  calculations  were  carried  out  for  54 

atmospheric  levels  and  two  sets  of  a  priori  statistics:  global  and  local  (with 

decreasing  all  covariance  matrix  elements  by  half).  The  examples  of  the 

calculation  results  are  given  in  Fig.4.1-4.10.  It  is  necessary  to  remind,  that 

when  calculating  the  error  matrices  (Pokrovsky  and  Timofeyev,  1972),  we  define  a 
potential  sensing  accuracy,  i.e.  minimal  attainable  error  for  examined  measurement 
error  levels.  Other  possible  factors,  which  can  increase  the  retrieval  errors 
(radiative  mode!  errors,  absolute  calibration  errors,  etc.)  are  not  taken  into 
account.  At  last,  let  us  note,  that  adduced  data  have  to  do  with  cloudless 
atmosphere  and  two  surface  types  -  land  and  sea.  Surface  emissivities  are  assumed 
to  be  known.  In  this  connection,  the  errors  analyzed  can  be  named  as  "potential" 
ones. 

Fig.4.1,  4.2  illustrate  the  potential  errors  of  vertical  temperature  profile 
retrieval  from  outgoing  radiance  measurements  in  1500  channels  and  from  SMM/T  and 
SMM/T-2  measurements,  respectively,  for  three  RMS  measurement  errors  in  the  terms 
of  brightness  temperature  (aT0=O.l;  0.3;  1.0K).  It  can  be  seen  from  the 
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Fig. 4.1. Potential  errors  of  the  temperature  profile 

retrieval  for  the  three  levels  of  measurement 


errors  ctTb. 

1500  measurements,  sea,  global  statistics. 
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Figures,  that  retrieval  errors  decrease  with  increasing  measurement  accuracy.  At 
0^=0. IK  the  complete  radiance  experiment  makes  possible  to  retrieve  T(p)  with  the 
errors  smaller  than  0.9K  in  the  wide  pressure  range  (10-1000  mb).  In  the  most  part 
of  this  pressure  range,  the  retrieval  errors  are  0.4-0.7K.  At  the  same  measurement 
error,  SSM/T+SSM/T-2  system  gives  the  maximal  retrieval  errors  equal  to  about  1.3K 
in  the  10-1000  mb  pressure  range. 

At  Ojb=1.0K,  the  complete  radiance  experiment  allows  to  reach  the  accuracies  equal 
to  0.5-1.4K  in  comparison  with  those  of  1.6-3.0K  for  SSM/T+SSM/T-2  system.  To 
compare  more  visually  the  retrieval  errors  of  two  experiment  designs,  those  are 
given  in  Fig.4.3  at  Ojb=0.3K  It  is  seen,  that  at  some  atmospheric  levels  (10-1000 
mb),  the  retrieval  advantages  of  the  complete  radiance  experiment  reach  to  1.0K 
and  to  larger  values  at  levels,  which  are  higher  than  10  mb.  The  last  is  a  result 
of  lacking  the  channels  with  high  optical  density  in  the  SSM/T+SSM/T-2  system.  The 
advantages  of  the  complete  radiance  experiment  for  the  10-1000  mb  layer  can  be 
clearly  seen  in  Fig.4.4. 

Thus,  the  calculations  of  potential  accuracy  of  atmospheric  temperature  sounding 

confirm  the  above  result,  that  outgoing  microwave  radiance  contains  the  noticeable 
additional  information  (as  compared  with  SSM/T+SSM/T-2  measurements)  about 

atmospheric  temperature  structure. 

An  influence  of  the  atmospheric  humidity  profile  variations  on  the  temperature 
sounding  accuracy  is  illustrated  by  Fig.4.5.  In  this  Figure,  two  profiles  of  T(p) 
retrieval  errors  are  given:  1  -  for  a  case  of  accounting  for  q(p)  variations;  2  - 
for  a  case  of  known  q(p).  The  effect  of  humidity  profile  variations  is  appreciable 
in  the  200-1000  mb  layer  and  can  reach  to  0.1-0.3K  for  the  complete  radiance 

experiment  and  to  about  0.5K  for  SSM/T+SSM/T-2  system  in  the  lower  troposphere. 

It  is  well-known,  that  accuracy  of  atmospheric  temperature  sounding  depends  on  the 
initial  guess  and  on  "width"  of  a  priori  statistics  used.  Fig.4.6  demonstrates  an 
influence  of  a  priori  statistics  on  the  accuracy  of  atmospheric  temperature 

sounding.  It  can  be  seen,  that  accuracy  of  the  atmospheric  temperature  retrieval 
by  use  of  "local"  statistics  at  aTg=0.3K  is  higher  by  0.2-0.4K  for  the  complete 
radiance  experiment  and  by  0.3-0.7K  for  SSM/T+SSM/T-2  system  .  The  smaller 
influence  of  a  priori  statistics  on  the  retrieval  with  using  1500  radiance 

measurements  is  caused  by  the  fact,  that  this  experiment  allows  to  obtain  a 

greater  information  about  the  atmospheric  temperature  state  and,  as  a  consequence, 
a  priori  statistics  slightly  effects  the  retrieval  results.  It  is  interesting  to 

note,  that  the  accuracy  of  temperature  profile  retrieval  at  the  10-400  mb 
atmospheric  layer  by  SSM/T+SSM/T-2  system  with  using  local  statistics  is  near  to 
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Fig. 4. 3. Potential  errors  of  the  temperature  profile 
retrieval  for  two  measurement  designs: 

1500  measurements  and  SSM/T+SSM/T-2. 

Sea,  global  statistics,  measurement  error— 0.3K. 
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Fig. 4. 4. Potential  errors  of  the  temperature  profile  retrieval 

in  the  10-1000mb  layer  for  two  njeasurement  designs: 
1500  measurements  and  SSM/T+SSM/T-2. 

Sea,  global  statistics,  measurement  error-0. 3K. 
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Fig.4. 5. Potential  errors  of  the  temperature  profile 
retrieval  for  the  fixed  and  variable  profiles 
of  the  atmospheric  humidity. 

Sea,  global  statistics,  measurement  error-0. 3K. 

1  -  1500  points,  variable  humidity 

2  -  1500  points,  fixed  humidity 

3  -  SSM/T+SSM/T-2,  variable  humidity 

4  -  SSM/T+SSM/T-2,  fixed  humidity 
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Fig.4. 6. Potential  errors  of  the  temperature  profile 
retrieval  using  global  or  local  statistics. 
1500  measurements  and  SSM/T+SSM/T-2. 
Sea,  measurement  error-0. 3K. 

1  -  1500  points,  global  statistics 

2  -  SSM/T+SSM/T-2,  global  statistics 

3  -  1500  points,  local  statistics 

4  -  SSM/T+SSM/T-2,  local  statistics 

i 


63 


that,  accessible  for  complete  radiance  experiment  using  a  global  statistics.  At 
the  lower  atmosphere,  the  complete  radiance  experiment  gives  the  errors,  which  are 
noticeably  smaller  than  those  for  SSM/T+SSM/T-2  measurements  even  for  global 
statistics. 

As  it  has  been  mentioned,  the  matrix  error  calculations  were  carried  out  for  54 
atmospheric  pressure  levels.  In  a  number  of  inverse  problem,  it  is  necessary  to 
know  the  mean  temperatures  of  some  atmospheric  layers.  Thus,  the  operative  system 
of  the  temperature  sounding  TOVS  performs  T(p)  retrieval  for  40  atmospheric 

pressure  levels  and  gives  the  information  about  the  mean  virtual  temperature 
(Reale  et  al.,  1994).  In  Fig.4.7,  we  reproduce  the  temperature  retrieval  errors 
for  two  measurement  designs  (at  aTQ=0.3K)  and  two  vertical  resolution.  The 

retrieval  accuracy  in  18  layers  is  higher  by  0.2-0.3K  for  both  experiment  designs 
in  comparison  with  the  attainable  for  54  atmospheric  levels.  Effect  of  vertical 
averaging  is  slightly  larger  for  the  complete  radiance  experiment.  It  is  evidently 
connected  with  the  larger  correlation  between  T(p)  retrieval  errors  at  different 
levels  for  SSM/T+SSM/T-2  system  having  only  12  channels  in  comparison  with  the 
complete  radiance  experiment  involving  1500  measurements.  It  is  important  to 

emphasize,  that  the  complete  radiance  experiment  makes  a  possibility  to  retrieve 
the  temperatures  in  18  layers  with  high  potential  accuracies  of  0.2-0.9K  at 

Otb=0.3K.  In  the  most  part  of  the  examined  10-1000  mb  atmospheric  layer,  the 
retrieval  errors  are  about  0.5K.  The  corresponding  errors  for  DMSP  sounding  system 
at  the  same  measurement  error  are  higher  by  0.5-1. OK. 

A  joint  influence  of  vertical  averaging  and  a  priori  statistics  on  the  accuracy  of 

the  mean  temperature  sounding  in  18  atmospheric  layers  at  OpB=0.3K  is  demonstrated 
in  Fig.4.8.  It  can  be  concluded,  that  complete  radiance  experiment  (1500 
measurements)  with  using  local  statistics  allows  to  achieve  the  very  high 

retrieval  accuracy  (0.2-0.6K)  for  average  temperatures  in  18  layers  of  the  10-1000 
mb  atmospheric  range.  The  same  retrieval  accuracies  for  DMSP  system  are  worse  by 
0.3-0.6K.  Above  this  atmospheric  range,  the  advantages  of  the  complete  radiance 

experiment  achieve  to  1-4K,  that  is  due  to  an  absence  of  optical  dense  measurement 

channels  in  the  SSM/T+SSM/T-2  system. 

An  influence  of  surface  emissivity  ^  on  the  temperature  retrieval  accuracy  is 

demonstrated  by  Fig.4.9.  This  impact  is  incidental  if  the  emissivity  is  known.  The 
exception  is  the  accuracy  of  temperature  retrieval  from  SSM/T+SSM/T-2  measurements 

near  underlying  surface.  In  this  case,  the  retrieval  accuracy  above  land  (e=0.95) 

is  significantly  better  than  that  above  sea.  It  is  connected,  in  our  opinion,  with 

special  shape  of  weighting  function  corresponding  to  the  "transparent"  measurement 
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Fig. 4. 7. Potential  errors  of  the  retrieval  of  the  temperature 
profile  on  54  levels  and  the  mean  temperature  in 
18  atmospheric  layers  for  two  measurement  designs: 
1500  measurements  and  SSM/T+SSM/T-2. 

Sea,  global  statistics,  measurement  error  —  0.3K. 

1  -  1500  points,  54  levels 

2  -  SSM/T+SSM/T-2,  54  levels 

3  —  1500  points,  18  layers 

4  -  SSM/T+SSM/T-2,  18  layers 
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Fig. 4. 8. Potential  errors  of  the  mean  temperature 
retrieval  in  18  atmospheric  layers  using 
global  or  local  statistics. 

1500  measurements  and  SSM/T+SSM/T-2. 

Sea,  measurement  error  -  0.3K. 

1  -  1500  points,  global  statistics 

2  -  SSM/T+SSM/T-2,  global  statistics 

3  -*•  1500  points,  local  statistics 

4  -  SSM/T+SSM/T-2,  local  statistics 
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channel  of  50.500  GHz  for  SSM/T  radiometer  above  sea  surface  (see  Fig.2.8). 

Fig.4.10  demonstrates  the  ratio  between  random  and  systematic  (null-space) 
components  in  T(p)  retrieval  overall  error  (see  Eqs.  (4.3)-(4.4))  for  examined 
experiment  designs  and  Otq=0.3K.  It  can  be  seen,  that  systematic  retrieval  error, 

caused  by  finite  thickness  of  the  temperature  weighting  functions,  gives  the 
maximal  contribution.  But  these  components  are  equal  to  0.3-1.0K  for  the  complete 
radiance  experiment,  and  to  0.8-1.8K  for  the  SSM/T+SSM/T-2  system. 

The  examples  of  similar  calculations  of  atmospheric  humidity  retrieval  accuracies 
for  different  designs  of  outgoing  microwave  radiance  measurements  are  given  in 
Fig.4.11-4.13. 

In  Fig.4.11  we  reproduce  the  curves  ^(p)  =  aq(p)/aq(p)  for  two  examined  experiment 
designs  at  oTb=0.3K.  (Here  aq(p)  and  aq(p)  are  a  priori  and  a  posteriori  RMS 
uncertainties  of  atmospheric  humidity  profiles).  By  definition,  W(p)  function 
characterizes  a  residual  uncertainty  in  retrieval  atmospheric  humidity  profile 

with  respect  to  a  priori  uncertainty.(  A  case  of  W(p)=0  corresponds  to  absolutely 
exact  retrieval  of  the  humidity  profile,  W(p)=l  -  to  the  lack  of  information  about 
q(p)  in  outgoing  radiance  measurements.).  Fig.4.11  demonstrates  the  essential 

advantages  of  the  complete  radiance  experiment  in  comparison  with  SSM/T+SSM/T-2 
system.  At  a  number  of  atmospheric  levels,  the  xIJ(p)  functions  corresponding  to 
examined  measurement  designs  differ  by  0.1-0.2.  If  we  turn  to  frequent-used 
relative  retrieval  error  aq(p)/q  (here  q  is  the  mean  atmospheric  humidity 
profile),  the  complete  radiance  experiment  and  SSM/T+SSM/T-2  system  give  the 
retrieval  errors  of  5-10%  and  10-30%,  respectively.  Cited  data  prove  the  above 
conclusion  that  SSM/T+SSM/T-2  measurements  does  not  make  possible  to  derive  all 
available  information  about  atmospheric  humidity  profiles  contained  in  outgoing 
microwave  radiance. 

Dependence  of  humidity  retrieval  errors  (in  the  terms  of  lF(p)  function)  on 

measurement  errors  is  given  in  Fig.4.12  for  the  complete  radiance  experiment.  The 
accuracy  of  humidity  sounding  increases  with  measurement  error  decrease.  It  is 
especially  observed  in  the  600-800  mb  layer,  that  is  connected  with  vertical 
behavior  of  humidity  weighting  functions. 

The  contribution  of  null-space  and  random  components  (at  oTb=0.3K)  into  overall 
error  of  atmospheric  humidity  profile  retrieval  (in  the  terms  of  W(p)  function) 
demonstrates  Fig.4.13.  As  well  as  for  the  case  of  atmospheric  temperature  profile 
retrieval  (see  Fig.4.10),  the  principal  contribution  is  introduced  by  null-space 
errors  caused  by  the  finite  thickness  of  humidity  weighting  functions. 
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Fig.4. 10. Potential  errors  (  overall,  null-space,  random  )  of 

the  temperature  profile  retrieval  for  two  measurement 
designs:  1500  measurements  and  SSM/T+SSM/T-2. 

Sea,  global  statistics,  measurement  error  —  0.3K. 

1  -  1500  points,  overall  error 

2  -  SSM/T+SSM/T-2,  overall  error 

3  -  1500  points,  null-space  error 

4  -  SSM/T+SSM/T-2,  null-space  error 

5  -  1500  points,  random  error 

6  -  SSM/T+SSM/T-2,  random  error 
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Fig. 4. 11. Function  ^  =  crq(p)/crq(p)  for  two  measurement 

designs:  1500  measurements  and  SSM/T+SSM/T-2. 
Sea,  global  statistics,  measurement  error  -  0.3K. 
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Fig. 4. 12. Function  'f/=<rq(p)/(rll(p)  at  different  measurement 
errors.  Sea,  global  statistics,  1500  measurements. 
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Fig. 4. 13. Potential  errors  (  overall,  null-space,  random  )  of 
the  humidity  profile  retrieval  for  1500 
measurements.  Sea,  global  statistics,  measurement 
error  -  0.3K. 
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4.3. Numerical  experiments  on  temperature-humidity  sounding 
of  the  atmosphere 

To  illustrate  the  potential  possibilities  of  the  atmospheric  temperature-humidity 
sounding  for  different  designs  of  outgoing  microwave  radiance  measurements,  we 
have  conducted  the  numerical  experiments  on  the  closed-loop  scheme.  The  complete 
cycle  of  the  numerical  experiments  includes  the  following  steps: 

1. For  given  model  of  atmospheric  physical  state,  i.e.  specific  vertical 
temperature  T(p)  and  humidity  q(p)  profiles,  the  brightness  temperatures  of 
outgoing  radiance  TB(f)  are  calculated  by  use  of  the  radiative  transfer  model, 
described  in  Chapter  2. 

2. The  measurement  errors  are  simulated  by  adding  the  E(f)  values,  produced  by 

random-number  generator  and  normally  distributed  with  assigned  variance  to 

the  TB(f)  values. 

3. The  joint  estimates  of  atmospheric  temperature  and  humidity  profiles  (T(p)  and 
q(p))  are  retrieved  by  iterative  method  (formula  (4.18)).  All  unknown  parameters 
of  the  surface,  namely,  emissivity  and  surface  temperature  are  derived 
simultaneously. 

4. The  retrieval  errors  are  determined  by  comparison  of  the  initial  atmospheric 
profiles  (T(p)  and  q(p))  with  the  estimates  (T(p)  and  q(p)). 

It  is  necessary  to  note,  that  when  solving  the  inverse  problem,  we  have  used  the 
above  iterative  method,  and  the  quasi-optimal  estimates  are  defined,  when  the 

iterative  process  converges. 

The  examples  of  the  temperature  retrieval  errors  (AT(p)=  T(p)-T(p))  for  different 
atmospheric  models  (ctTb=0.3K;  the  retrievals  are  carried  out  on  54  levels)  are 
given  in  Fig.4. 14-4.17  for  two  measurement  designs  (the  complete  radiance 

experiment  and  the  SSM/T+SSM/T-2  system).  The  advantages  of  the  complete  radiance 
experiment  in  the  temperature  sounding  are  obvious  on  the  most  atmospheric  levels. 
The  .errors  of  the  temperature  profile  retrieval  in  the  10-1000  mb  atmospheric 
layer  do  not  exceed,  as  a  rule,  IK  for  the  complete  radiance  experiment.  Those  for 
the  SSM/T+SSM/T-2  system  may  reach  to  2-3K. 

The  influence  of  outgoing  radiance  measurement  designs  on  the  accuracy  of  humidity 
profile  retrieval  is  illustrated  by  Figures  4.18-4.21.  These  Figures  represent  the 
initial  q(p)  profiles,  the  estimates  q(p)  retrieved  from  the  measurements 
simulating  the  complete  radiance  and  DMSP  experiments,  and  the  mean  atmospheric 
profile  q(p).  The  advantages  of  the  complete  radiance  experiment  are  especially 
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Fig. 4. 14. Temperature  retrieval  errors  for  two  measurement 
designs:  1500  measurements  and  SSM/T+SSM/T-2 
(numerical  experiment). 

Land,  measurement  error  -  0.3K. 
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Fig. 4. 15. Temperature  retrieval  errors  for  two 

measurement  designs:  1500  measurements 
and  SSM/T+SSM/T-2  (numerical  experiment). 
Land,  measurement  error  -  0.3K. 

1  -  1500  measurements 

2  -  SSM/T+SSM/T-2 
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Fig.  16. Temperature  retrieval  errors  for  two 

measurement  designs:  1500  measurements 
and  SSM/T+SSM/T-2  (numerical  experiment). 
Land,  measurement  error  -  0.3K. 

1  -  1500  measurements 

2  -  SSM/T+SSM/T-2 
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Fig. 4. 17. Temperature  retrieval  errors  for  two 

measurement  designs:  1500  measurements 
and  SSM/T+SSM/T-2  (numerical  experiment) 
Sea,  measurement  error  -  0.3K. 

1  -  1500  measurements 

2  -  SSM/T+SSM/T-2 


300 


Fig. 4. 18. Humidity  profile  retrievals  for  two  measurement 
designs:  1500  measurements  and  SSM/T+SSM/T-2 
(numerical  experiment). 

-Land,  measurement  error  -  0.3K. 
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Fig. 4. 19. Humidity  profile  retrieval  for  two 

measurement  designs:  1500  measurements, 
and  SSM/T+SSM/T-2  (numerical  experiment). 
Land,  measurement  error  -  0.3K. 

1  -  the  mean  profile 

2  -  exact  profile 

3  -  1500  measurements 

4  -  SSM/T+SSM/T-2 
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Fig. 4. 20. Humidity  profile  retrieval  for  two  measurement 

designs:  1500  measurements  and  SSM/T+SSM/T-2. 
Sea,  measurement  error  -  0.3K. 

1  -  the  mean  profile 

2  -  exact  profile 

3  -  1500  measurements 

4  -  SSM/T+SSM/T-2 
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Fig. 4. 21. Humidity  profile  retrieval  for  two  measurement 

designs:  1500  measurements  and  SMM/T+SMM/T-2 
Land,  measurement  error  -  0.3K. 

1  -  the  main  profile 

2  -  exact  profile 

3  -  1500  measurements 

4  -  SSM/T+SSM/T-2 


revealed  in  Fig.4.19,  4.20.  It  is  seen,  that  the  complete  radiance  experiment 
allows  to  retrieve  a  relatively  sophisticated  vertical  structure  of  q(z)  profiles 
with  high  accuracy.  The  remote  sensing  system  of  DMSP  satellite  gives  strongly 
smoothed  q(p)  and,  as  a  consequence,  the  errors  of  the  humidity  retrieval  from  the 
DMSP  data  are  essentially  larger. 

Similar  examples  of  the  numerical  experiments  on  atmospheric  temperature 
retrievals  in  18  atmospheric  layers  are  given  in  Figures  4.22-4.24.  The  advantages 
of  the  complete  radiance  experiment  are  evident.  In  the  most  cases,  the 
atmospheric  temperature  profile  may  be  retrieved  using  the  complete  radiance 
experiment  data  with  the  errors  which  are  less  than  IK. 
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Fig. 4. 22. Errors  of  the  temperature  retrievals  in  18  atmospheric 

layers  for  two  measurement  designs:  1500  measurements 
and  SSM/T+SSM/T-2  (numerical  experiment). 

Land,  measurement  error  -  0.3K. 

1  -  1500  measurements 

2  -  SSM/T+SSM/T-2 
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Fig. 4. 23. Errors  of  the  temperature  retrievals  in  18  atmospheric 
layers  for  two  measurement  designs:  1500  measurements 
and  SSM/T+SSM/T-2  (numerical  experiment). 

Land,  measurement  error  -  0.3K. 

1  —  1500  measurements 

2  -  SSM/T+SSM/T-2 
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Fig. 4. 24. Errors  of  the  temperature  retrievals  in  18  atmospheric 
layers  for  two  measurement  designs:  1500  measurements 
and  SSM/T+SSM/T-2  (numerical  experiment). 

Sea,  measurement  error  -  0.3K.  ‘  ■ 

1  —  1500  measurements 

2  -  SSM/T+SSM/T-2 
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5.PRINCIPAL  RESULTS  AND  CONCLUSIONS 

The  studies  performed  in  the  frame  of  the  Special  Contract  SPC-94-4086  have  given 
the  following  results: 

1. A  new  code  (MIT-SPbSU  program)  for  calculations  of  outgoing  microwave  brightness 
temperatures  TB(f)  for  a  case  of  cloudless  atmosphere  has  been  developed.  Liebe’s 
(1993)  algorithm  of  accounting  for  atmospheric  radiative  characteristics  is 
assumed  as  a  basis  of  this  program.  Distinctive  characteristic  of  this  program  is 
a  possibility  to  account  the  variational  derivatives  of  brightness  temperature 
with  respect  to  different  physical  parameters  of  atmosphere  and  underlying  surface 
(atmospheric  vertical  temperature  and  humidity  profiles,  surface  emissivity). 

Performed  comparisons  of  TB  values  calculated  by  MIT-SPbSU  code  and  independent 
data  have  demonstrated  a  high  quality  of  the  developed  program.  Dependencies  of  TB 
variational  derivatives  with  respect  to  different  parameters  of  atmospheric  state 
have  been  analyzed.  It  is  shown,  for  example,  that  a  problem  of  atmospheric 
humidity  remote  sounding  is  essentially  nonlinear. 

2. The  software  meant  for  the  solution  of  the  inverse  problems  and  for  the 
calculations  of  error  estimates  of  temperature-humidity  satellite  sounding  in  the 
1-200  GHz  spectral  range  has  been  developed.  This  code  is  based  on  the  known 
iterative  method  (Polyakov  and  Rozanov,  1989;  Marks  and  Rodgers,  1993),  but  it  has 
been  modified  to  retrieve  the  optimal,  in  the  statistical  sense,  estimates  of  the 
solution. 

3. The  calculations  of  covariance  matrices  K  of  outgoing  microwave  radiance  for 

tbtb 

limited,  but  representative,  ensemble  of  atmospheric  states  corresponding  to 
global  observations  have  been  conducted.  1500  measurement  channels  in  the  1-200 
GHz  spectral  range  with  special  channel  choice  are  used  to  describe  in  detail  the 
optical  atmospheric  characteristics  and  to  simulate  "the  complete  radiance 
experiment".  The  calculations  of  K  are  carried  out  for  two  types  of  underlying 

TBTB 

surface  -  sea  and  land.  The  choice  of  the  emissivity  models  employed  for 
atmospheric  studies  has  been  justified  by  emissivity  calculations  for  available 
part  of  examined  range  using  the  special  developed  code. 

4. The  eigenvalues  and  eigenvectors  of  outgoing  radiance  covariance  matrices  have 
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been  calculated  for  two  experiment  designs  -  "complete  radiance  experiment"  and 
experiment  using  SSM/T  and  SSM/T-2  measurements  from  DMSP  satellite.  According  to 
Kozlov’s  criterion  (Kozlov,  1966),  a  number  of  independent  information  components 
contained  in  outgoing  radiance  measurements  has  been  defined  for  different 
Measurement  designs  and  error  levels.  It  is  shown,  that  for  examined  atmospheric 
parameter  set,  corresponding  to  the  global  observations,  a  number  of  independent 
measurement  components  ranges  from  15  to  31  depending,  first  of  all,  on  random 
measurement  errors  Op  and  surface  type.  So,  for  example,  at  oTn=0.5K  there  is  a 
possibility  to  "pick  out"  18  (land)  and  22  (sea)  independent  information 
components,  caused  by  temperature  and  humidity  variations,  from  outgoing  microwave 
radiance  spectra.  These  values  essentially  exceed  the  number  of  SSM/T  and  SSM/T-2 
channels,  intended  for  temperature  and  humidity  sounding. 

The  additional  calculations  have  demonstrated  that  the  most  part  of  independent 
information  components  is  connected  with  temperature  profile  variations.  But,  the 
humidity  profile  variations  form  6-8  independent  information  components  in  the 
space  of  measured  signals. 

Analysis  of  the  eigenvalues  of  the  radiance  covariance  matrices,  corresponding  to 
SSM/T  and  SSM/T-2  radiometer  measurements  shows,  that  there  are  2-3  "redundant" 
measurement  channels  in  SSM/T  radiometer  at  aTe=0.5-1.0K,  and  the  same  one  channel 
for  SSM/T-2  radiometer  at  ctXb=1.0K.  This  fact  (with  regard  for  performed 
information  content  analysis  for  "complete  radiance  experiment")  justifies 

indirectly  to  non-optimal  positions  of  DMSP  spectral  channels. 

5.Information  content  of  outgoing  microwave  radiance  measurements  has  been  defined 
by  two  criterions:  according  to  Shannon  and  Kozlov.  These  calculations 

demonstrated  the  "information  advantage"  of  the  complete  radiance  experiment  in 
comparison  with  SSM/T  and  SSM/T-2  measurements  even  with  accounting  for  different 
number  of  measurement  channels. 

6.0n  .the  basis  of  the  error  matrix  calculations,  the  temperature  and  humidity 
retrieval  errors  have  been  calculated  for  different  experiment  designs, 
measurement  errors,  "global"  and  "local"  statistics,  and  54  levels  and  18 

atmospheric  layers.  Analysis  of  the  results  has  confirmed  the  principal 

conclusions  received  on  the  basis  of  information  content  calculations: 
a).Accuracy  of  the  temperature  profile  retrieval  from  the  complete  radiance 

experiment  data  is  essentially  better  than  that  for  SSM/T  and  SSM/T-2 

measurements.  In  particular,  at  aTe=0.3K,  the  retrieval  errors  for  experiment  with 
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1500  channels  are  smaller  by  1.0K  for  a  case  of  using  global  statistics. 

b)  .Temperature  retrieval  errors  decrease  noticeably  for  both  experimental  designs 
if  average  temperatures  in  18  atmospheric  layers  are  defined  and/or  the  local 
atmospheric  statistics  is  used  for  the  retrieval.  In  this  case,  for  example,  at 
c*pB=0.3K,  the  complete  radiance  experiment  makes  a  possible  to  decrease  a 
"potential"  error  to  0.3-0.9K  in  the  10-1000  mb  layer.  The  minimal  temperature 
retrieval  errors  which  can  be  achieved  using  DMSP  sounding  system  are  greater  by 
0.7-1.0K. 

c) .Similar  advantages  of  the  complete  radiance  experiment  are  observed  during 
humidity  atmospheric  sounding.  At  ctTb=0.3K,  minimal  potential  errors  of  relative 
humidity  profile  retrieval  are  5-15%  and  10-30%  in  the  300-1000  mb  atmospheric 
layer  for  1500  measurements  and  SSM/T  +  SSM/T-2  measurements,  respectively. 

7.The  numerical  simulation  of  the  satellite  experiments  on  joint  atmospheric 
temperature  and  humidity  profile  retrieval  has  been  carried  out.  The  analysis  of 
the  retrieval  errors  confirms  the  principal  conclusions  obtained  as  a  result  of 
the  information  content  and  potential  accuracy  studies. 

Thus,  performed  numerical  studies  of  information  content  of  outgoing  microwave 
measurements  and  the  accuracy  of  atmospheric  temperature  and  humidity  sounding  let 
draw  the  following  main  conclusions: 

l.Outgoing  microwave  radiance  spectra  in  the  1-200  GHz  range  contain  a 
considerable  quantity  of  information  about  vertical  temperature  and  humidity 
atmospheric  profiles.  In  the  terms  of  independent  information  components,  this 
information  is  characterized  by  15-31  parameters  depending  on  the  levels  of 
experimental  errors.  The  main  part  of  these  components  describes  the  temperature 
variations  (11-21  parameters),  the  other  (5-11  parameters)  characterize  the 
humidity  profile  variations. 

2.A  number  of  measurement  channels  of  SSM/T  +  SSM/T-2  system  of  DMSP  satellite  is 
deficient  to  extract  all  existing  information  about  atmospheric  temperature  and 
humidity  profiles  from  outgoing  microwave  radiance  spectra.  Besides,  the  channel 
spectral  locations  are  not  optimal  in  a  number  of  cases. 

So  the  most  important  problem  is  to  compute  the  optimal  experiment  designs  for 
microwave  radiometer  with  about  30  channels  on  the  basis  of  recent  methods  of 
measurement  condition  optimization. 
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